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(N ' ABSTRACT 

We present a method to analyze the spectral energy distributions (SEDs) of 
O ' young stellar objects (YSOs). Our approach is to fit data with pre-computed 

2-D radiation transfer models spanning a large region of parameter space. This 
CN ' allows us to determine not only a single set of physical parameter values but 



the entire range of values consistent with the multi-wavelength observations of a 



given source. In this way we hope to avoid any over-interpretation when modeling 
^ • a set of data. We have constructed spectral energy distributions from optical 

to sub-mm wavelengths, including new Spitzer IRAC and MIPS photometry, 
for 30 young and spatially resolved sources in the Taurus-Auriga star-forming 
region. We demonstrate fitting model SEDs to these sources, and find that we 
correctly identify the evolutionary stage and physical parameters found from 
previous independent studies, such as disk mass, disk accretion rate, and stellar 
temperature. We also explore how fluxes at various wavelengths help to constrain 
physical parameters, and show examples of degeneracies that can occur when 
fitting SEDs. A web-based version of this tool is available to the communitjQ. 
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stars: formation — stars: pre-main-sequence 
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Introduction 



Studying the circumstellar environment of YSOs is essential to our understanding of 
the pre-main-sequence evolution of stars. Since it is not possible to observe a single YSO 
through various stages of evolution over thousands or millions of years, this has to be 
done statistically, by observing large numbers of YSOs, and inferring from these possi- 



decade (e.g. 


Tamura e 


t al. 


1991 




Kenvon et al. 


1993 




Burrows et al. 


1996; 


Whitnev et al. 


1997 




Close et al. 


1997 




Lucas & Rochel 


1997: I 


3taDelfeldt et al. 


1998: 


Padeett et al. 


1999: 


Cotera et al. 


2001 


), giving us some insight into the evolution of low-mass YSOs in nearby 



star-forming regions. However, only a limited number of YSOs can be directly resolved 
(~ 100), most of which only marginally. To obtain a quantitative picture of both low and 
high-mass star-formation, we require observations of a larger sample, by studying the tens 
of thousands of unresolved YSOs seen in nearby star-forming regions (e.g Taurus-Auriga, 
Perseus, Orion) as well as more distant star-forming regions (e.g. M16, M17, NGC6334). 

In order to study unresolved or close-to unresolved YSOs in these distant regions, we can 
resort to mult i- wavelength photometry, from which we construct SEDs. The main question 
we seek to answer in this series of papers is: how can we make the most of the information 
contained in the SEDs while limiting any over-interpretation? 

Many methods have been sugg ested and u sed in the past to classify an d interpret YSO 



SEDs, such as spectral indi ces (e.g. lLadalll9870 or color-color diagrams (e.g. lLada fc Adams 



19921 : lAUen et al.ll2004l ). In iRobitaille et al.l (120061 . hereafter Paper I), we presented a grid of 
200,000 model SEDs spanning a large range of possible evolutionary stages and stellar masses, 
and using these models, we explored what could be learned from using spectral indices and 
color-color plots. The main advantage of these two techniques is that they make it easy and 
fast to classify the SEDs of sources, although it is not necessarily trivial to extract information 
relating to the physical conditions in these objects. Furthermore, many YSOs now have data 
at more than four wavelengths, and although color-color or color-magnitude plots with more 
than two dimensions can be constructed to make use of this extra information, only sources 
that have data in all of the required filters can be used. In large-scale mult i- wavelength 
studies of star-formation, it is not uncommon to have a significant number of sources that 
lack at least one datapoint. Therefore, it is usually not possible to make the best use of all 
of the data available for each source using these techniques. 

Another approach to analyzing SEDs is to compute radiation transfer models, first 
assuming a given circumstellar dust and gas geometry, as well as dust properties, predicting 
the emergent SED, and finding a set of parameters that reproduce the observations - it is a 
typical inverse problem. The main advantage of this technique is that one gets an insight into 
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the actual physical properties of the source, rather than simply quantifying the shape of the 
SED. Another advantage of using SED fitting is that it can make use of any data available, 
while not being limited by the lack of a datapoint at a given wavelength. Of course, the 
more data are present, the better the parameters will be constrained. 

If a source has been observed at many different wavelengths from optical to mm wave- 
lengths, it is probable, though not certain, that the set of parameters for the model SED 
that will match the data well will be fairly unique, provided that all the parameters actually 
affect the shape of the SED, and that there are no degeneracies. However, many star- forming 
regions have been covered by several surveys, providing flux measur ements at only half a 



dozen to a dozen wavelengths (e.g . 2MASS, Smtzerl RAC & MIPS; ISkrutskie et all 12006 



Fazio et al. 2004; Rieke et al. 2004 



Werner et al.ll2004 ) for large numbers of sources. When 



studying a large number of YSOs, one wants to know not only one set of physical parameters 
that can explain the data for each source, but all the different sets of parameters that can 
explain the data for each source, in order to avoid any over-interpretation. Fitting SEDs to 
these many sources by trial and error would be problematic: it is very likely that there is 
not a unique set of parameters that can explain the data available for a given source, and it 
would be extremely time-consuming to explore parameter space manually for each source, 
let alone thousands of sources. 

Our proposed solution to this problem is to precompute a large set of radiation transfer 
models that provide a reasonably large coverage of parameter space. We can then use this 
set of models to find all the combinations of parameters that can explain the data for a given 
source. 

This paper describes the method that we use to fit these models to data for single sources, 
and presents examples of what can be learned from using such a technique. In Section [2] 
we describe the technical details of the fitting method. In Section [3] we demonstrate the 
use of the fitting method on Taurus- Auriga sources as a proof-of-concept: we first describe 
the sample of sources and the data used to construct the SEDs ( §3.11) ; then we compare the 
values of a selection of physical parameters with independent estimates ( §3.21) ; and finally 
we show how different combinations of datapoints constrain physical parameters for a given 
source, and discuss degeneracies that can arise ( §3.3]) . Additionally, we show an example of 
a resolved source for which we can measure fluxes in apertures smaller than the source itself, 
and demonstrate that we can fit the SEDs from the different apertures simultaneously ( §3.4[) . 
An example of large-scale study of a star-forming region will be presented in Robitaille et 
al. (2007, in preparation) and in Indebetouw et al. (2007, in preparation). 
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2. Technical overview 

The SED fitting tool that we have developed uses the 200,000 YSO model SEDs (20,000 
sets of physical parameters and 10 viewing angles) presented in Paper I. The models consist 
of pre-main-sequence stars with different combinations of axisymmetric circumstellar disks, 
infalling flattened envelopes, and outflow cavities, covering a wide range of parameter space. 
The main caveat for this work is that we assume that stars of all masses form via accretion 
through a disk and envelope. However, we note that the grid of models covers a large range 
of parameter space, so that we assume as little as possible about the specifics of the accretion 
physics. 

The parameter ranges covered by the models span those determined from observations 
and theories. The parameters we vary fall into three categories: the central source parame- 
ters (stellar mass, radius and temperature), the infalling envelope parameters (the envelope 
accretion rate, outer radius, inner radius, cavity opening angle and cavity density), and the 
disk parameters (disk mass, accretion rate, outer radius, inner radius, flaring power, and 
scaleheight). The stellar masses M* are sampled between 0.1 and 50 M & , the stellar ages 
£* are sampled between 10 3 and 10 7 yr, and the s tellar radii R± and temperatures T+ are 



found from M* and t* using evolutionary tracks (IBernasconi fc Maederl 1 19961 ; ISiess et al. 



20001 ) . The disk and envelope parameters are then sampled randomly within ranges that 
depend on the age of the central source. For example, the disk mass was sampled from 
Mfak/Mi, ~ 0.001 — 0.1 at early evolutionary stages, and a wider range of masses extending 
down to Mdi s k/M* = 10~ 8 between 1 and lOMyr to allow for the disk dispersal stage. The 
dust grain models used in the radiation transfer models are the following: the densest regions 
of the disk (nn 2 > 10 10 cm~ 3 ) use a grain model with a size distribution that decays expo- 
nentially for sizes larger than 50 fim extending up to 1 mm, and the rest of the circumstellar 
geometry uses a grain size distribution with an average particle size slightly larger than the 
diffuse ISM, and a ratio of total-to-selective extinction Ry = 3.6. We assume a gas-to-dust 
ratio of 100 (note that the results can be scaled to different gas-to-dust ratios since only the 
dust is taken into account in the radiation transfer). For more details about the models and 
the ranges of parameters sampled, including caveats, we refer readers to Paper I. 

We have convolved these models with common filter bandpasses ranging from optical 
(e.g. UBVRI), to near and mid-IR (e.g. 2MASS JHK, Spitzer IRAC & MIPS), far-IR, and 
sub-mm (e.g. IRAS, SCUBA). The details of the convolution of the model SEDs with the 
filters used in this paper is described in Appendix |A] A large range of filters is available in 
the online fitting tool, and additional filters can be added on request from users. 

Let us assume that we have photometry of a source at N wavelengths Aj (j=i...N), and 
that the flux densities are F u (Xi) with uncertainties cr(Aj). In addition, we assume that the 
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fluxes were measured in apertures A(\i). Finally, we assume that the source lies in a distance 
range <i m i n to c? max . 

In order to fit convolved model fluxes to data we first interpolate the fluxes to the 
apertures A(Xi) used to perform the photometry, for a number of distances d between d min 
and G? max , and scale them to the appropriate distance. We are able to do this since each SED 
was computed for 50 apertures between 100 and 100,000 AU. 

At each of these distances, we fit the convolved model fluxes to the data using optimal 
scaling, allowing the visual extinction Ay to be a free parameter. Writing the convolved, 
interpolated, and scaled model fluxes as M„(A;), the overall pattern P„(A*) that is being fit 
to the data is then 



P„(Ai) = M„(Ai) x i -°- 4A v^Mv 



X 



where K\ i and K y are given by an extinct i on law . We use an extinction model calculated 
with the method of iKim. Martin, fc Hendry! (119941 ) that fits a typical Gala c tic IS M curve 
modified for the mid-IR extinction properties derived by llndebetouw et al. J2005h (Wolff, 
private communication). Since P u (\) is nonlinear, it is convenient to work with log [F u ] and 
log [P v ] rather than F v and P v so as to transform this into a linear problem. The model that 
is fit to the data is then: 



log 10 [P v (\i)] = log 10 [M„(\i)] - 0AA 



Ky 



V" 



(2) 



where the free parameter Ay can be determined from optimal scaling. The unbiased 
fluxes and flux variance in log 10 space are given by: 



>gio \ F A) 
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(3) 
(4) 
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Once the free parameter Ay has been determined, we compute the x 2 -P er_ datapoint 
value of the fit: 
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(logio [Fu(Xi)] ) - logio [PJM 
a((log 10 [FJ(Ai)» 



) 



2 



(6) 



where iV is the number of datapoints being fit. Throughout this paper, the x 2 values 
mentioned are per datapoint. 

This fitting process is repeated for a range of distances between <i min and <i max , and a 
best-fit d, Ay, and x 2 value can be found for each model SED. The process is repeated for 
all the models in the grid. In this way, all the model SEDs in the grid are compared with 
the data, and the parameters of the model SEDs that fit a source well can be analyzed. 

The process described above does not make any assumptions about the spatial extent 
of the source. For example the apertures may be smaller than the extent of the source (see 
§3.41 for an example). Usually however, one measures the total flux of a source by using an 
aperture larger than the source. With the knowledge that the source is not larger than a 
given aperture, we can cut down the number of model SEDs that fit its observed SED well. 

For example, if a source in Taurus is unresolved at a given wavelength Aj, then we 
should not be fitting models that would have been resolved at that particular wavelength. 
More generally, if we know that the apparent extent of a source is always smaller than the 
apertures that the fluxes were measured in, then we want to eliminate all models whose 
apparent extent is larger than the aperture in at least one wavelength. Note that we do not 
want to eliminate all models that are physically larger than the aperture, but only those 
that appear larger than the aperture (the apparent size can be smaller than the physical 
size). To use this information, we compute the wavelength-dependent outermost radius 
i?Ai(S = So/2) at which the surface brightness S of the model is equal to half of the peak 
surface brightness E . We can then reject all models that have -Ra^S — ^o/2) > A(X i ) in 
at least one wavelength Aj. 

The online version of the fitting tool is hosted on a dedicated Web serveJE At the 
moment, it is possible to fit only a single source at a time, but we plan to upgrade this to 
allow users to upload data-files in order to fit multiple sources in late-2007. 



2 http :/ /www. astro .wise . edu/protostars 
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3. Proof-of-concept using known Taurus- Auriga sources 



3.1. The data 



In this section we analyze whether by fitting model SEDs to observed SEDs we are 
able to correctly identify the evolutionary stage of YSOs, as well as determine the value of 
individual physical parameters. To do this we have fit the SEDs of well studied and resolved 
YSOs in the Taurus-Auriga star-forming region. Over the last few decades, Taurus has been 
the best studied star-forming region, and is in that respect the region for which the most 
data are available. We hav e constructed a sample of 30 sources by selecting those from 
Kenyon fc Hartmannl (119951 ) that have been spatially resolved as of mid-2006. The latter 
requirement ensures that we have a good a-priori knowledge of the evolutionary stage of the 
objects from direct observations (e.g. young protostars with infalling envelopes, or disks). 
FS Tau and FS Tau B were not included due to their small angular separation, which makes 
it difficult to construct separate SEDs, and UZ Tau E was not included because of insufficient 
SED coverage. DG Tau and DG Tau B were included, despite their small angular separation, 
as IRAC and MIPS data were available, allowing us to construct two separate well defined 
SEDs. 

A few sources are known fr om direct imaging t o be more comple x than axisymmetric 
objects (e.g. IRAS 04016+2610: IPadeett et"aDll999l : Iwood et alJboOlh . but it is interesting 
to see what we can derive from fitting the SEDs of such sources. Indeed, when looking at 
more distant star-forming regions, we have no a-priori knowledge of the complexity of the 
geometry of an unresolved source, and it is interesting to know what we can learn using 2-D 
models 

We have compiled the SED s for these 30 sources using optical (UBVRI), near- (JHKL) 
and mid-IR (MNQ) data from iKenvon fc Hartmannl (119951 ). J HK fluxes from t he 2M ASS 



point-source catalog (ISkrutskie et al.ll2006l ). IRAC data from lHartmann et al 



2005 



1) and 



Weaver fc Jones 



Luhman et al.l (120061 ). IRAS data from the IRAS point-source catalog and from 
(119921 ) (the latter was used over the former when available), SHAR C II 350 /xm, SCUBA 
450um, and 850 fim sub-mm data from the compilation presented in Andrews fc Williams 
J2005h . and CSO 624 and 769 /xm data from beckwith fc Sargentl jl99lh . 



In ad dition, we have measu red the IRAC fluxes for the sources not presented in lHartmann et al. 



(120051 ) or iLuhman et al.l (120061 ). and re-measured the IRAC fluxes of RY Tau and DG Tau, 
which are saturated in all bands. The IRAC data (PI Fazio - P00006 & P00037) and the 
MIPS data (PI D. Padgett - P03584) were retrieved from the Spitzer Science Center Archive. 
To measure the fluxes of sources saturated in IRAC, we used the publicly available iracworks 
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codej written by Tom Jarrett. Finally, we have performed PSF photometry of MIPS 24 and 
70 /im data, using a custom photometry code. For sources that were mildly saturated in 
MIPS, we ignored the affected pixels and ensured that the outer parts of the PSF were being 
correctly fit. 

The list of sources along with the data are shown in Tables [BJ [U [TO], and [11] in Ap- 
pendix IB] 

When fitting the SEDs, it is technically possible to use all the data available for each 
source, but in cases where several measurements were available at the same wavelength, we 
used the highest quality one. For example, if MIPS data are available, we used it instead of 
IRAS, since IRAS has a poorer resolution, and could result in confusion. Similarly, space- 
based IRAC data are usually preferred over ground based mid-IR data. We applied the 
following rules: 

• when 2MASS JHK data was available, it was used instead of previous JHK measure- 
ments. 

• when IRAC 3.6 /im data was available, the L-band flux was not used. 

• when IRAC 4.5 and/or 5.8 /im data were available, the M-band flux was not used. 

• when MIPS 24 /im data was available, the IRAS 25 /im flux was not used. 

• when MIPS 70 /im data was available, the IRAS 60 /im and 100 /im data were not 
used. 

• in general, flux measurements were always preferred over upper limits. 

Additionally, we decided not to use the N and Q data when fitting the SEDs, due to large 
uncertainties in the fluxes, filter profiles, calibration, and zero-magnitude fluxes, but we quote 
these flux values in the tables for reference. A lower limit of 25% was imposed on the flux 
uncertainties for optical, L & M, IRAS, and sub-mm data, so that any flux measurement 
with a smaller uncertainty saw its uncertainty increased. This was done to account for 
variability and uncertainties in the zero-magnitude fluxes in the optical, uncertainties in the 
filter profiles and zero magnitude fluxes in the L & M bands, and uncertainties in the absolute 
flux calibration for IRAS and sub-mm wavelengths. A lower limit of 40% was imposed on the 
flux uncertainties for the CSO sub-mm data, as we did not have the transmission profiles for 



3 ht t p : / / spider . ipac . caltech . edu /staff/j arrett /irac 
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these observations: since the full width at half maximum (FWHM) of the filters used for the 
769 yum observations is large (190 /xm), differences in the filter profile would be noticeable. 
For the transmission profiles of these two bands we used gaussians centered at 624 and 
769 /im with FWHM 67 and 190 /xm respe c tively , which we convolved with the atmospheric 
transmission profile used in iDowell et al.l (120031 ). Additionally, a lower limit of 10% was 
imposed on the remaining datapoints, to account mainly for uncertainties in the absolute 
calibration and photometry extraction. 



3.2. Comparison of derived physical properties to other methods 

The only assumptions we made when fitting the observed SEDs with our model SEDs 
was that the sources were all situated in a distance range of 120 — > 160AU (to rule out too 
luminous or too faint sources) , that the foreground interstellar extinction was no more than 



Ay = 20 (IWhittet et al.ll200ll ). and we used the condition that none of the sources appeared 
larger than the apertures used to measure the fluxes. The aperture radii assumed for the 
data are listed in Table [12] in Appendix [Bl 

In Figure [1] we show the SEDs for all the sources along with the best fitting model for 
each source, and all the models with y 2 — xlest < 3 (where Xmm i s the x 2 P er datapoint of the 
best fitting model for each source). Although this cutoff is arbitrary, it provides a range of 
acceptable fits to the eye. A purely statistical argument would of course be more desirable, 
but the sampling of our models in 14 dimensional space is too sparse to resolve well enough 
the shape of the minima in the chisquared 'surface' to derive formal confidence intervals. 
Using a 'chi-by-eye' approach does carry the risk of overestimating the uncertainties on the 
various parameters, but we argue that this is better than under-estimating the uncertainties, 
which would lead to over-interpretation. Furthermore, intrinsic variability of the sources, 
along with the fact that no young stellar object will be as perfectly symmetrical and well 
behaved as our models, means that any 'formal' confidence interval would likely be too strict. 

All sources except DM Tau and GM Aur have their SED reasonably well fit by our 



model SEDs. These sources are known to show a near-IR deficit in their SEDs ( Rice et al. 



20031 ; ICalvet et al.l 120051 ) . suggesting that the inner region of their disks has been cleared of 
dust. Although we do include a substantial number of models with inner holes in our grid, 
our models systematically over-estimate the observed mid-IR fluxes. This may be because 
we evacuate the disks completely below a given radius, rather than having small but non-zero 
amounts of dust inside the 'hole'. We plan to address this in a future grid of models. This 
is an example where fitting model SEDs to data provides us with feedback to improve our 
models. 
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The SED for CoKu Tau/1 is fit by only one model within the goodness-of-fit constraints. 
The model that does fit well is that of an edge-on YSO with a remnant accreting envelope 
(see Tabled]). Based on the SED fit alone, an edge-on inclinat ion seems re a sonab le, as the 



SED is double-peaked, a typical signature of an edge-on disk. IStark et all (120061 ) modeled 
spatially resolved near-IR images of CoKu Tau/1, and found that it does indeed requires 
a low-mass envelope, suggesting that CoKu Tau/1 is more evolved than a typical Class I 
source. However, they find that the images are best-fit using a 64° inclination rather than 
edge-on. This seems more plausible as the central source is visible in these images (the 
central source would not be visible in an edge-on disk). The central source is a binary 
system with separation 0.24" (~33AU at 140 pc), which will likely have evacuated some 
material from the inner disk. At this point, we can only speculate about the reasons for 
the discrepant results between the SED and image fitting: Perhaps one of the sources in the 
binary is surrounded by a small close to edge-on disk inside the larger circumbinary disk. 
This may explain why one of the sources in the binary is redder in NICMOS images than 
the other source. We plan to further model the multi- wavelength images of CoKu Tau/1 to 
test this hypothesis. 

In general, we note that the main effect a young binary system will have on its cir- 
cumstellar material will be to clear out, at least partially, the innermost regions of the disk 
or envelope of gas and dust. In fact, this is the primary reason that large inner holes are 
included in our models. The main issue will likely be the determination of the properties of 
the central sources. In cases where the binary separation is large enough (e.g. CoKu Tau/1 
or GG Tau), the geometry will become more complex, with a circumbinary disk, and possibly 
small individual disks around each source. In these cases, the disk parameters will also be 
wrong, as the geometry is not being modeled correctly. In the future we plan to explore in 
more detail the effects of a binary system on the determination of the physical parameters 
of YSOs. 

In the remainder of this section, we compare a selection of physical properties of the 
sources derived from fitting model SEDs to the sources with values quoted in the literature: 



Evolutionary Stage and Viewing Angle The first result that we can compare with 
known values is whether the correct evolutionary stage is identified, i.e. whether young 
protostars with infalling envelopes or disks are correctly identified as such (see Table [TJ). 
To do this, we look at the range of envelope accretion rates that provide a good fit to each 
source, and specifically the lower value of the range of good-fitting models. If this lower value 
is greater than zero, then this tells us that the source cannot be explained by a disk-only 
source within the goodness-of-fit constraints applied. Conversely, if the lower value is zero, 
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then this tells us that the source can be explained by a disk with no infalling envelope. In 
most cases where the lower value is zero, the upper value is non-zero, which means that it is 
impossible to rule out that there may be very optically thin envelopes surrounding the disks. 
We also list the ranges of inclinations to the line-of-sight that provide a good fit as there is in 
some cases a degeneracy between evolutionary stage and inclination. We can make several 
remarks about the results: 

• All known disk sources (marked as 'Disk' in Table [T]) except DM Tau are well fit by 
models with no infalling envelope. Furthermore, the maximum value of the accretion 
rate providing a good fit for these sources does not exceed lO~ 6 M yr -1 . This is 
expected, as we showed in Paper I that envelopes with accretion rates lower than 
~ 10 -6 M Q yr _1 do not contribute significantly to the SED. 

• Most sources that are known to have infalling envelopes (marked as 'Embedded' in 
Tabled]), cannot be fit by disk-only models. 

• The exceptions to the above point are DG Tau B, IRAS 04016+2610 and IRAS 04248+2612, 
which are also known to have infalling envelopes, but can be fit by disk-only models as 
well as models with infalling envelopes. However, for all three sources, the disk-only 
models that provide good fits are all viewed edge-on (the lower limit on the range 

of inclinations for disk-only models, i.e. models for which M env = 0, which fit these 
sources well is larger than 80°). This is a good demonstration that even using an 
SED with data ranging from optical or near-IR to sub-mm is not always enough to 
unambiguously distinguish between various geometries once viewing angle is taken into 
account. 

In practice, one could eliminate model fits based on statistical arguments. For example, 
in a star-forming cluster, observing a source with an edge-on disk has a low probability both 
because one is less likely to observe a disk with an edge-on inclination, and because such a 
source would be much fainter. As discussed further in this section, the edge-on disk models 
that fit the observed SEDs well require the luminosity of the central star to be larger than 
100 Lq, which is unrealistic, as we would in this case expect to observe many more non 
edge-on sources with such luminosities. 

We note that in most cases, the inclination of the line of sight is not well determined. 
The only constraints we find are that none of the disk-only sources are seen exactly edge-on 
(87° in Table 1), and that CoKu Tau/1 cannot be fit by any non edge-on models as discussed 
previously. 
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Stellar Temperature In Table [2] and in Figure [2] we compare the temperatures derived 
from SED fitting to those corresponding to the spectral types listed in iKenyon fc Hartmann 
( 119951 ). The latter were determined from spectroscopic observations, and therefore provide 
an independent and more accurate measure of the temperature. We find that for all sources 
but one, the best-fit value is close to the known value. The exception to this is UY Aur, but 
we note that the range of temperatures for this source (3540 — > 8090 K) is consistent with 
the literature value (4060 K). 

AB Aur is the only source in our sample that is known to have a high temperature 
(~ 10, 500 K), and we note that it is also the source for which the best-fit model has the 
highest temperature (11,767K). Generally, we find that we have correctly identified the 
temperatures of the sources to better than ±0.2 orders of magnitude. 

Disk Mass We c o mpare the disk masse s der ived from SED fitting to y alues listed in 



Dutrey et al.l (119961 ) , iKitamura et al.l (120021 ) , and lAndrews k, Williams! (120051 ) . These values 
were determined from sub-mm and mm observations. In Table [3] we list the various literature 
values, including an average value, a minimum value defined to be the smallest value of 
Mdisk — o"(^disk) quoted, and a maximum value defined to be the largest value of M disk + 
cr(M disk ) quoted. We compare this range to the range found from SED fitting (this is also 
shown graphically in Figure [2]). Our values assume a gas-to-dust ratio of 100. As can be 
seen from Figure [2], for the very young sources (marked 'Embedded' in Tabled]) we do not 
constrain the disk mass well. For example IRAS 04361+2547 could, on the basis of the 
SED, have a disk mass between 2.5 x 10~ 5 M Q and 3.2 x 10~ 2 M Q . This is because, in the 
early stages of evolution, when the disk is deeply embedded inside the infalling envelope, 
the relative contributions of the disk and envelope to the SED are difficult to disentangle. 
For disk-only sources (marked 'Disk' in Table [1]), we find that our values generally agree 
with the literature values to within less than an order of magnitude. This is fairly accurate 
considering that the disk masses in our model grid range over nine orders of magnitude, 
and that literature values often rely on simple models (e.g. power laws), and assume that 
the sub-mm flux arises from an optically thin isothermal region of disks. Furthermore, the 
choice of a dust opacity law will also affect the disk mass obtained. 



Disk Accretion Rates In Table 0] and in 
derived from SED fitt ing to values listed in 



Hartmann et al.l (119981 ). and iMohanty et al.l (120051 ). Most of these literature values were 



igure [2] we c o mpare the disk accretion rate s 



Valenti et al.l (Il993l ). lHartigan et al.l (Il995[ ) 



derived from UV and optical spectroscopy, and therefore represent an independent and more 
accurate estimate of the disk accretion rate than SED fitting. We find that the agreement 
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between values derived from SED fitting and literature values is reasonable over the three 
orders of magnitude spanned. 

The values obtained from SED fitting appear to be slightly larger than those taken 
from the literature. This could be due to a small inconsistency in our models where the 
accretion luminosity inside the dust destruction radius but outside the magnetic truncation 
radius is, for simplicity, being emitted as a stellar photon, with a resulting stellar spectrum. 
Thus, to match a given near-IR excess, a slightly higher accretion rate is required. For 
examp le, RY Tau is the s ource for which our estimate deviates the most from literature 



values; lAkeson et al.l (120051 ) found that for this source, the accretion luminosity from the gas 
disk inside the dust destruction radius contributed significantly to the total disk accretion 
luminosity. Future versions of our grid of models will include disk emission inside the dust 
destruction radius. 



Disk/Envelope Inner Radius All the sources in the sample can be fit by disks and en- 
velopes with no inner holes (-R^in = -^sub where R su b is the dust sublimation radius defined in 
Paper I), with the exception of DM Tau, GM Aur, IRAS 04302+2247, and IRAS04325+2402 
(see Table [5]). As described in Paper I, the inner radius is the same for the disk and the 
envelope. For most sources, the upper limit to the disk inner radius is larger than i? su b, and 
it is generally not possible to rule out the existence of holes in any of the sources. However, 
in the case of DM Tau, GM Aur, IRAS 04302+2247, and IRAS04325+2402, the lower limit 
being larger than R sn b does rule out disks or envelopes with no inner holes. GM Aur has 
previously been found to require an inner hole of order ~ 4AU in order to explain its SED 



( IRice et al.ll2003l ). More recently, Spitzer IRS spectroscopy has confirmed that both GM Aur 



and D M Tau show a near -IR deficit, suggesting inner hole sizes of 24 AU and 3AU respec- 



tively (jCalvet et al.ll2005l ). We fit the SED of GM Aur with models having an inner hole 



size between ~ 1 and 24 AU, and DM Tau with models having an in ner hole size between 



3.4 and 20 AU: these values are consistent with the values found by lCalvet et al. 



The above comparisons show that in most cases, the parameters of the best fit model 
agree well with previously published values. Even in cases where the best-fit value does not 
agree with the published values, the range of values providing a good fit is consistent with 
previously published values. 

In addition to the parameters presented above, we also present the ranges of values 
providing a good fit for the stellar mass, bolometric luminosity, disk outer radius, and disk 
scaleheight h at 100 AU. The values are presented in Tables [5] and [6j The disk outer radii 
are generally fairly uncertain, and range from tens to hundreds of AU. The stellar masses 
are likely to be constrained because we determine temperatures well, and our model grid 
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uses evolutionary tracks, relating the stellar temperatures to the stellar masses. We do not 
have literature values to compare the derived scale heights to. However we do expect the 
derived values to be consistent with the true values: as shown in Paper I, the disk flaring 
power does have an effect on the SED, meaning that we should be somewhat sensitive to the 
disk scaleheight. The values are poorly constrained for embedded sources, where the disk 
does not contribute significantly to the SED. Finally, the main uncertainty in the bolometric 
luminosity is due to uncertainties in the inclination. For example the sources for which the 
luminosity is very poorly constrained (from a few L Q to over 100 L ), such as DG Tau B 
or IRAS 04016+2610, are those for which we cannot distinguish between an edge-on disk 
or a young protostar with an infalling envelope. These two possibilities would result in two 
very different bolometric luminosities. If the inclination can be constrained, then we can 
potentially determine the true luminosity of the source very accurately from SED fitting. 

3.3. Constraining parameters: uniqueness of fits, degeneracies, and 

wavelength coverage 

As mentioned in Section [H when fitting model SEDs to multi-wavelength observations 
of a source, it is important to have a grasp not only of which set of parameters provides 
a good fit, but also whether it is the only set of parameters that does so. For example, a 
source may be well fit by a model with a 1 M Q central source and a 300 AU disk with a 
mass of Mdisk = 0.01 M & disk, but it may be that the radius of the disk does not actually 
have an effect on the SED at the wavelengths of the available data, and that any disk radius 
would provide a suitable fit. Furthermore, two completely different sets of parameters may 
fit data equally well (e.g. the sources from Section 13.21 that can be fit by SED models of 
young protostellar objects or older edge-on disks). 

The reader may wonder how we can claim to determine 14 parameters from fitting SEDs 
to data that may or may not span a large wavelength range. The answer is that we do not 
claim to determine any of the parameter values, but we aim to find how well constrained 
each parameter is. For example, from fitting IRAC data alone, we may find that none of 
the parameters are well constrained and that adding an extra datapoint only allows us to 
determine for example an upper limit on one of the parameters. 

In this section, we show how the YSO model SED fitting tool can be used to analyze how 
well various parameters are constrained when using different combinations of datapoints. In 
addition, we show examples of degeneracies that can arise in parameter space. This is by no 
means an exhaustive study, and we encourage readers to use the online fitting tool to explore 
in further detail how well various parameters are constrained for a specific set of data. 
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The first example that we show is AA Tau, which is a known T-Tauri source with optical 
variability. Figure [3] shows the model SED fits to the data, and a selection of parameters for 
these fits (disk mass and accretion rate, envelope accretion rate, and stellar temperature). 
For each combination of datapoints we show all models with \ 2 ~ Xmin < 3 as in Section 
13.21 As before, we assume a distance range of 120 — > 160pc. We first fit the SED using 
only IRAC points, and we find that 11,674 model SEDs satisfy the goodness of fit criterion 
(~ 5% of all model SEDs!). The disk mass, disk accretion rate and envelope accretion rate 
are not at all constrained apart from very high values of the temperature being ruled out 
(this is because models with very high temperatures would be too luminous to explain the 
IRAC fluxes). This suggests that using IRAC fluxes alone does not necessarily provide a 
good estimate of the evolutionary stage of an object. Adding JHK data reduces the number 
of good fits to 2,826 model SEDs, but does not provide any significant improvement in the 
determination of the four parameters that we show. Adding IRAS 12 /xm, MIPS 24 /zm and 
MIPS 70 jum does have an important effect which is to provide a much more strict upper 
limit on the envelope accretion rate (~ 10 il^yr -1 ), and also provides a lower limit on 
the disk mass of 1O _6 M . Furthermore, the disk accretion rate and the central source 
temperature are better determined. This echoes our findings presented in Paper I that data 
beyond ~ 20 /im is very helpful in constraining the evolutionary stage of a source. Adding 
the optical data rules out models that have both low central source temperatures and low 
disk accretion rates, leaving a degeneracy between the optical data being explained by a 
higher temperature or a higher disk accretion rate. Finally, adding the sub-mm data has the 
predictable effect of providing a strong constraint on the disk mass, placing it at ~ 1O _2 M . 
The degeneracy between central source temperature and disk accretion rate is resolved, and 
the disk accretion rate is estimated at ~ 10 -8 — lO _7 M /yr. 

Our second example is IRAS 04361+2547, which is a known protostar surrounded by an 
infal ling envelope and with molecular outflo ws, as shown by spatially resolved observations 
(e.g. lTerebey et al.lll990uTamura et al.lll99ll ). Figure H] shows the model SED fits to the data, 
and a selection of parameters for these fits (disk mass, envelope accretion rate, and stellar 
mass and temperature). As for AA Tau we show all models with \ 2 ~ Xmm < 3. We first fit 
the SED using only IRAC points, and we find that 1,959 model SEDs satisfy the goodness of 
fit criterion. As before, the disk mass and envelope accretion rate are not at all constrained. 
Stellar mass and temperature are also very poorly constrained. Adding JHK data reduces 
the number of good fits to 183, and does provide an improvement in the determination of the 
envelope accretion rate. However, we note that there is a bimodal distribution of models in 
(M env , M disk ), with one peak centered at fairly high accretion rates (~ lO~ 5 M /yr), and one 
centered at zero accretion rates (which are mostly edge-on disk-only models). A bimodal 
distribution is also seen in (T*,M*). The models with high temperatures correspond to 
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the edge-on disk-only models, which require on average a larger central source luminosity 
to match the observed fluxes, while the models with lower temperatures correspond to the 
young models with large envelope infall rates. This is a good example of a young protostar - 
edge-on disk degeneracy. The addition of IRAS 12 /zm and IRAS 25 /mi reduces the number 
of edge-on disk models that fit well, but does not rule them out. Adding IRAS 60 fim and 
IRAS 100 /im rules out all edge-on disk models, and in the process rules out all models with 
M* > 1M . Finally, adding the sub-mm data further constrains the parameters further, 
albeit not much. By this stage, the SED fit is fairly well constrained. The envelope accretion 
rate is constrained to within an order of magnitude, the disk mass is determined to within 
two orders of magnitude, the temperature is determined to be no more than 4,000 K, and 
the stellar mass is determined to be no more than 1 M & . 



3.4. Analysis of resolved sources 

As mentioned in Section [2J when analyzing the SED of a resolved object, it is possible 
to specify apertures smaller than the apparent extent of the source. In this section, we 
demonstrate this by modeling IRAS 04368+2557 using our model SED fitter and using 
solely IRAC data. 

In order to model IRAS 04368+2557, we found the total flux in the four bands in six 
different apertures: 4" , 8" , 16" , 32" , 48" , and 64" . A composite color image of the source is 
shown in Figure [5j along with the position of the six apertures. 

We imposed a distance constraint of 120 to 160pc to the object, as well as an Av to the 
object of to 20. We did not impose any conditions such as an inclination of the object to 
the line of sight, nor the cavity size, as we wish to find whether the fitting tool can determine 
such parameters solely from the aperture information. 

The left plots in Figure M shows the fluxes for the six apertures along with the best 
fitting aperture-dependent model SED. Also shown are the near-IR, far-IR, and sub-mm 
datapoints, but these are initially not used in the fitting process. Using the parameters 
of the best fit, we produced a model image, shown in Figure [7J by rerunning this model 
with a higher signal-to-noise. The model SED fitter correctly finds an SED that matches 
the data adequately in all apertures, determines the type of object (infalling envelope), the 
inclination to the line of sight (edge-on), and the colors of the object. The cavity angle is 
not a perfect match, but this is not surprising since we only provided the fitting tool with 
circularly averaged fluxes, containing no information about any possible asymmetries. 

We note however that all the models that fit the resolved IRAC data well systematically 
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overestimate fluxes at far-IR wavelengths (as seen in the top left plot of Figure [6]). The plots 
on the right of Figure [6] show that fitting the data at other wavelengths simultaneously 
deteriorates the fit at IRAC wavelengths. Furthermore, our models do not reproduce the 
peak i n luminosity at the center of the object (as seen in Figure [7]). It has recently been 



found (lLoinard et al.ll2002l ) that the central source of IRAS 04368+2557 is a binary star with 
a projected separation of 25 AU, meaning that we expect a large cavity at the center of the 
infalling envelope. As mentioned previously, although our model grid does include models 
with inner holes, we evacuate the holes completely, rather than leaving low levels of dust. 
A small but non-zero amount of dust inside an inner hole may be enough to provide the 
IRAC fluxes observed, while decreasing the amount of far-IR flux. The parameters of the 
best-fitting models for the IRAC data only, and for the full SED, are listed in Tabled 

As well as demonstrating SED fitting of resolved sources, this is an illustration of how 
detailed modeling of an individual source can be used to improve our models for a future 
grid. In this future grid we also plan to produce images at a number of different wavelengths 
for all the models in the grid, allowing us to fit multi-wavelength images of resolved YSOs 
without first averaging the flux in circular apertures. 



4. Summary and Conclusions 

We have developed a method to fit observed YSO SEDs using a pre-computed grid of 
models. Although we vary the values of 14 physical parameters in our models (c.f. Paper I), 
we do not claim to be able to determine all of these parameters from fitting observed SEDs. 
Instead, we are interested in determining which parameters can be contrained, if at all, and if 
so, what range of values provide acceptable fits. We have compiled data from UBVRI to sub- 
mm wavelengths for 30 YSOs in the Taurus-Auriga star-forming region. The evolutionary 
stages, stellar temperatures, disk masses, and disk accretion rates derived from fitting the 
model SEDs to the data are in good agreement with independently determined values (for 
example from spectroscopy). In cases where the best fitting model is not in agreement with 
the literature value, the range of parameter values of the model SEDs that provide a good 
fit is in general still consistent with this value. We have also demonstrated how adding 
fluxes at various wavelengths helps constrain different parameters. For example we find, 
as in Paper I, that data in the range 20-100 fim in addition to shorter wavelength data is 
very useful in constraining parameters such as the envelope accretion rate, and thus the 
evolutionary stage of a source. An online fitting tool making use of this method is available 
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to the community^. In future, we plan to fit spectra (e.g. Spitzer IRS data) and polarization 
measurements simultaneously with the broadband SED, in order to place further constraints 
on the geometry and the chemistry of the circumstellar environment of YSOs. 
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Fig. 1. — SEDs for the 30 Taurus- Auriga sources analyzed in this paper. Filled circles are the 
flux values listed in Tables [HJ [9J [10] and [TT1 Triangles are upper limits. Error bars are shown if 
larger than the data points. The solid black line indicates the best fitting model, and the gray lines 
show all models that also fit the data well (defined by \ 2 ~ XjLst < 3, where \ 2 1S the value per 
datapoint). The dashed line shows the SED of the stellar photosphere in the best fitting model. 
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Fig. 1. — Continued. 
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Fig. 1. — Continued. 
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Fig. 2. — Comparison of the values of stellar temperature (top left), disk accretion rate (top 
right), and disk mass (bottom left and right) found from SED fitting with values quoted 
in the literature and derived from different methods. The dashed line shows where the 
values would be equal, and the gray area shows where the agreement is better than ±0.2 
orders of magnitude for the temperature, and ±1 order of magnitude for the disk mass and 
accretion rate. Error bars in the y-direction indicate the range of values of the models shown 
in Figure [TJ Error bars in the x-direction show the range between the lowest value and 
the highest value quoted in the literature. For the stellar temperature, only one reference 
( iKenyon &: Hartmannlll995l ) was used, and therefore no uncertainties are shown on the x- 
axis. The likely uncertainty quoted by the authors is ±1 spectral class, corresponding roughly 
to ±100K. 
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Fig. 3. Left: model SED fits to the observed SED of AA Tau using from top to bottom: IRAC points only, JHK + IRAC 

points, JHK + IRAC + MIPS 24 //m + MIPS 70 /urn, UBVRI + JHK + IRAC + MIPS 24 /jm + MIPS 70 ftm, and UBVRI + 
JHK + IRAC + MIPS 24 jim + MIPS 70 fim + sub-mm points (when more than 200 model SEDs fit the data well, we show 
the 200 best fit SEDs, and one in ten SEDs beyond this). Center and right: a selection of parameters for the model fits. Filled 
circles arc the good fits shown in the SED plots. The grayscale shows the parameter space of the grid of models (smoothed). 
Each shade of gray shows an increase in a factor of 10 in the density of models. 
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Fig. 4. Left: model SED fits to the observed SED of IRAS 04361+2547 using from top to bottom: IRAC points only, 

JHK + IRAC points, JHK + IRAC + IRAS 12 & 25 /im, JHK + IRAC + All IRAS bands, and JHK + IRAC + All IRAS 
bands + sub-mm points (when more than 200 model SEDs fit the data well, we show the 200 best fit SEDs, and one in ten 
SEDs beyond this). Center and right: a selection of parameters for the model fits. Filled circles are the good fits shown in the 
SED plots. The grayscale shows the parameter space of the grid of models as for Figure [3] The alignment of all the fits along 
a curve in the T* vs plots (right) is due to the sampling of the models using evolutionary tracks, i.e. for a given M* and t* 
we find T* using evolutionary tracks (see Paper I). 
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Fig. 5.— Left: Three color image of IRAS 04368+2557 using IRAC 3.6 /an (blue - to 
SMJy.sr" 1 ), IRAC 4.5 /mi (green - to eMJy.sr" 1 ), and IRAC 7.8 /mi (red - to 12MJy.sr- 1 ). 
Right: The same image with the six photometry apertures overplotted (4", 8", 16", 32", 48", 
and 64"). 
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Fig. 6. — Top left: best-fit model to the IRAC data in the six apertures for IRAS 04368+2557 
(only the IRAC data is used). The IRAC datapoints are shown as filled black circles. The 
other datapoints are shown as gray circles to indicate that these were not used in the fitting. 
The apertures assumed for these points are those listed in Table [12] as before, i.e. 100" for 
JHK, 10" for MIPS 24/im, 60" for IRAS 12 & 25 /im, 120" for IRAS 60 & 100 /im, and 30" 
for the sub-mm data. For the IRAC datapoints, the fluxes in six apertures are shown (4", 
8", 16", 32", 48", and 64"), with the faintest fluxes corresponding to the smallest apertures, 
and the brightest fluxes corresponding to the largest apertures. The model SED in the 
different apertures is shown as open circles. The different colors correspond to the different 
apertures. Top right: best-fit model using the near-IR, far-IR, and sub-mm data in addition 
to the IRAC data in the six apertures. The bottom left and right panels show the same SED 
fits in the IRAC range. 
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Fig. 7.— Image of the best-fit model to the IRAC data only for IRAS 04368+2557, shown 
in the same bands (IRAC 3.6, 4.5, and 7.8 /im) and with the same scaling as the image in 
Figure [51 
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Table 1. Comparison of values of the envelope accretion rate found from SED fitting with 
the evolutionary stage of the object determined from resolved observations. 









SED Fitting 


values 






- U J 


i ( M 
I {IVlenv 


-> U ) 


DUUICt; IN exilic 


n/vui. oidgc 


A/Tin 




Best 




Max 




1V11I1 


IVlclX 


A/Tin 
IV 11 11 


iviax 


AA Tau 


Disk 












7.69x10- 


-8 


18 


81 


18 


63 


AB Aur 


Hae 












1.41x10" 


-6 


32 


63 


32 


57 


BP Tau 


Disk 







5.18x10 


-8 


7.69x10- 


-8 


18 


81 


18 


81 


CI Tau 


Disk 












5.72x10- 


-8 


18 


76 


18 


63 


CoKu Tau 1 


Embedded/Disk 


6.88x10 


-8 


6.88x10 


-8 


6.88x10- 


-8 






87 


87 


CY Tau 


Disk 

















18 


81 






DG Tau 


Disk 







5.48x10 


-9 


1.11x10" 


-7 


18 


81 


49 


81 


DG Tau B 


Embedded 







1.00x10 


-6 


1.00x10" 


-6 


87 


87 


18 


87 


DL Tau 


Disk 

















18 


63 






DM Tau 


Disk 


7.08x10 


-9 


6.19x10 


-7 


8.89x10- 


-6 






18 


49 


DN Tau 


Disk 












5.18x10" 


-8 


18 


81 


18 


63 


DO Tau 


Disk 







7.22x10 


-9 


7.22x10- 


-9 


32 


41 


18 


81 


DR Tau 


Disk 







7.22x10 


-9 


5.67x10- 


-7 


18 


76 


18 


76 


FT Tau 


Disk 







4.73x10 


-8 


7.21x10" 


-6 


18 


81 


18 


76 


GG Tau 


Disk 












2.31x10" 


-6 


18 


57 


18 


41 


GM Aur 


Disk 












1.01x10- 


-6 


18 


57 


18 


41 


HL Tau 


Embedded 


1.87x10 


-6 


o 1 o v 1 n 
Z.loX 1U 


-6 


1 £0 vin- 
i.oy x iu 


-5 






18 


32 


IQ Tau 


Disk 












5.18x10" 


-8 


18 


81 


18 


81 


IRAS 04016+2610 


Embedded 







9.28x10 


-7 


4.86x10- 


-6 


81 


87 


18 


41 


IRAS 04169+2702 


Embedded 


3.70x10 


-7 


1.52x10 


-6 


2.78x10- 


-5 






32 


76 


IRAS 04248+2612 


Embedded 







1.73x10 


-6 


5.10x10" 


-6 


81 


81 


18 


87 


IRAS 04302+2247 


Embedded 


1.65x10 


-6 


1.23x10 


-5 


1.23x10" 


-5 






18 


81 


IRAS 04325+2402 


Embedded 


1.20x10 


-6 


1.20x10 


-6 


1.51x10" 


-5 






32 


76 


IRAS 04361+2547 


Embedded 


8.51x10 


-7 


1.46x10 


-6 


3.00x10- 


-5 






18 


81 


IRAS 04365+2535 


Embedded 


7.16x10 


-7 


9.28x10 


-7 


2.64x10- 


-5 






18 


76 


IRAS 04368+2557 


Embedded 


9.11x10 


-6 


2.82x10 


-5 


5.05x10- 


-5 






41 


81 


L1551 IRS5 


Embedded 


5.54x10 


-6 


6.47x10 


-5 


2.96x10- 


-4 






18 


76 


LkCa 15 


Disk 












3.76x10- 


-7 


18 


70 


18 


63 


RY Tau 


Disk 












9.52x10- 


-8 


18 


76 


57 


76 


UY Aur 


Disk 












5.88x10- 


-7 


70 


70 


57 


81 



Note. — All accretion rates are M Q /yr. 'Embedded' refers to sources that are still surrounded by an infalling 
envelope, and 'Disk' refers to sources that are surrounded only by a circumstellar disk. The last four columns 
show the ranges of inclinations i for the disk-only models (M env = 0) which provide a good fit as well as the range 
of inclinations i for the embedded models (M env > 0) which provide a good fit. In the Monte-Carlo radiation 
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transfer code (see Paper I), photons are binned into ten viewing angles. The values quoted here correspond to 
the central value of the bins. The central values of the bins are, from face-on to edge-on respectively: 18°, 32°, 
41°, 49°, 57°, 63°, 70°, 76°, 81°, and 87° (to the nearest degree) 
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Table 2. Comparison of values of the central source temperatures from the literature with 

those found from SED fitting. 







SED Fitting values 


Source Name 


T 1 

A star 


T ■ 

± min 


Tbest 


T 

± max 


AA Tau 


4060 


3060 


4458 


4859 


AB Aur 


10500 


4881 


11767 


13452 


BP Tau 


4060 


3138 


3427 


5415 


CI Tau 


4060 


3486 


4255 


4859 


CoKu Tau 1 




4826 


4826 


4826 


CY Tau 


3720 


3855 


4329 


4329 


DG Tau 


4350^5080 


4314 


4549 


10722 


DG Tau B 




2580 


2706 


13618 


DL Tau 


4060 


4255 


4255 


4458 


DM Tau 


3720 


2762 


2993 


3045 


DN Tau 


3850 


3427 


4314 


4458 


DO Tau 


3850 


4123 


4123 


5209 


DR Tau 


4060 


3935 


4123 


8061 


FT Tau 




2783 


3060 


5013 


GG Tau 


4060 


3320 


4430 


4859 


GM Aur 


4730 


3303 


5126 


5126 


HL Tau 


4060 


2561 


o a aa 
3409 


A AO A 

4030 


IQ Tau 


3785 


3427 


4255 


5612 


IRAS 04016+2610 




2552 


2847 


13587 


IRAS 04169+2702 




2585 


2718 


4152 


IRAS 04248+2612 


3580 


2627 


2909 


4956 


IRAS 04302+2247 




2762 


2986 


3241 


IRAS 04325+2402 




2586 


2741 


3616 


IRAS 04361+2547 




2585 


3225 


4064 


IRAS 04365+2535 




2637 


2847 


4228 


IRAS 04368+2557 




2932 


3855 


3869 


L1551 IRS5 


4060^6030 


3070 


3651 


4900 


LkCa 15 


4350 


3017 


4150 


4204 


RY Tau 


5080 


4604 


4604 


5713 


UY Aur 


4060 


3540 


8090 


8090 



Note. — All temperatures are in K. 

References. - All literature spectral types were taken 
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from iKenvon fc Hartmann ( 1995), except t hat for DG Tau 



which is from IWhite fc Hillenbrand (|2004l ). We use Table 
A5 in IKenvon fc Hartmann! (|1995T ) to convert the spectral 



types to temperatures. 



-36- 



Table 3. Comparison of values of the disk mass from the literature with those found from 

SED fitting. 









Literature values 






SED Fitting 


values 


bourcG JN ctme 




rvUZ 


A05 


Mean 


IV /Tin 


l\/To-v- 


Min 


Best 


A/To v 


AA Tau 


-1.8 


-1.6 ± 0.3 


-1.9 ± 0.1 


-1.7 


-2.0 


-1.3 


-2.1 


-1.7 


-1.5 


AB Aur 






-2.4 ± 0.1 


-2.4 


-2.4 


-2.3 


-1.5 


-1.5 


-1.0 


BP Tau 


-2.9 


. . . 


-1.8 ± 0.1 


-2.0 


-2.9 


-1.7 


-3.0 


-1.9 


-1.4 


CI Tau 


-1.2 




-1.6 ± 0.1 


-1.4 


-1.7 


-1.2 


-2.0 


-1.5 


-1.4 


CoKu Tau 1 






~ -2.7 


-2.7 


-2.7 


-2.7 


-3.6 


-3.6 


-3.6 


CY Tau 


-1.8 


-1.3 ± 0.4 


- -2.2 


-1.5 


-2.2 


-0.8 


-2.1 


-1.7 


-1.7 


DG Tau 


-1.6 




-1.6 ± 0.1 


-1.6 


-1.7 


-1.6 


-3.7 


-1.2 


-1.2 


DG Tau B 




. . . 










-3.7 


-2.6 


-1.1 


DL Tau 






-1.1 ± 0.1 


-1.0 


-1.2 


-1.0 


-1.7 


-1.5 


-1.5 


DM Tau 


-1.6 


-1.7 ± 0.3 


-1.6 ± 0.1 


-1.6 


-2.0 


-1.4 


-3.1 


-2.4 


-2.2 


DN Tau 


-1.6 


-9.0 ± 2.5 


-1.5 ± 0.1 


-1.6 




-0.8 


-2.1 


-1.5 


-1.5 


DO Tau 


-1.7 


-2.7 ± 0.5 


-2.2 ± 0.1 


-2.0 




-1.7 


-1.9 


-1.5 


-1.5 


DR Tau 




-2.3 ± 0.3 


-1.7 ± 0.1 


-1.9 


-2.7 


-1.7 


-2.0 


-1.5 


-1.2 


FT Tau 


-1.9 




-1.9 ± 0.1 


-1.9 


-1.9 


-1.8 


-3.0 


-2.1 


-1.5 


GG Tau 






-0.7 ± 0.2 


-0.6 


-0.8 


-0.5 


-1.9 


-1.5 


-1.5 


GM Aur 


-1.4 


-1.5 ± 0.3 


-1.6 ± 0.1 


-1.5 


-1.9 


-1.2 


-2.6 


-1.5 


-1.4 


HL Tau 




-1.3 ± 0.1 


-1.2 ± 0.1 


-1.2 


-1.4 


-1.1 


-1.9 


-1.5 


-1.1 


10 Tau 




-1.4 ± 0.2 


-1.7 ± 0.1 


-1.5 


-1.7 


-1.2 


-2.3 


-1.5 


-1.4 


IRAS 04016+2610 






- -1.7 


-1.7 


-1.7 


-1.7 


-2.2 


-1.8 


-0.9 


IRAS 04169+2702 






- -1.5 


-1.5 


-1.5 


-1.5 


-3.0 


-1.8 


-1.3 


IRAS 04248+2612 






- -2.3 


-2.3 


-2.3 


-2.3 


-4.3 


-2.2 


-1.4 


IRAS 04302+2247 






- -1.5 


-1.5 


-1.5 


-1.5 


-5.5 


-2.7 


-2.2 


IRAS 04325+2402 






- -2.1 


-2.1 


-2.1 


-2.1 


-4.6 


-2.8 


-2.5 


IRAS 04361+2547 






- -2.0 


-2.0 


-2.0 


-2.0 


-4.6 


-3.6 


-1.5 


IRAS 04365+2535 






- -1.5 


-1.5 


-1.5 


-1.5 


-3.4 


-1.8 


-1.2 


IRAS 04368+2557 






~ -1.4 


-1.4 


-1.4 


-1.4 


-4.2 


-3.6 


-1.8 


L1551 IRS5 






- -0.3 


-0.3 


-0.3 


-0.3 


-3.5 


-1.5 


-0.9 


LkCa 15 




-1.4 ± 0.4 


-1.3 ± 0.1 


-1.3 


-2.1 


-1.0 


-2.1 


-2.0 


-1.7 


RY Tau 




-3.4 ± 1.0 


-1.8 ± 0.1 


-1.9 




-1.7 


-2.4 


-1.5 


-1.0 


UY Aur 






-2.8 ± 0.1 


-2.7 


-2.8 


-2.7 


-2.9 


-1.4 


-1.4 



Note. — All disk masses are in log [M Cl 



References. 



D96: iDutrev et all |l996l ). K02: iKitamura et~all |2002h 



A05: 



Andrews fc Williams! |2005l 
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Table 4. Comparison of values of the disk accretion rate from the literature with those 

found from SED fitting. 



Literature values SED Fitting values 



Source Name 




V93 


H95 


H98 


M05 


Mean 


Min 


Max 


Min 


Best 


Max 


AA Tau 




-8.15 


-6.90 


-8.48 


-8.56 


-8.02 


-8.56 


-6.90 


-7.96 


-7.16 


-6.8 


AB Aur 


















-7.13 


-6.92 


-5.88 


BP Tau 




-7.61 


-6.80 


-7.54 


-7.71 


-7.42 


-7.71 


-6.80 


-8.26 


-7.19 


-6.71 


CI Tau 




-7.83 


-6.80 


-7.19 




-7.27 


-7.83 


-6.80 


-7.31 


-6.93 


-6.58 


CoKu Tau 1 


















-10.08 


-10.08 


-10.08 


CY Tau 






-8.20 


-8.12 


-8.16 


-8.16 


-8.20 


-8.12 


-8.52 


-6.95 


-6.95 


DG Tau 






-5.70 




-6.93 


-6.32 


-6.93 


-5.70 


-11.2 


-6.07 


-6.07 


DG Tau B 


















-9.98 


-6.98 


-5.64 


DL Tau 




-7.63 


-6.70 






-7.17 


-7.63 


-6.70 


-7.16 


-6.93 


-6.93 


DM Tau 




-8.54 




-7.95 




-8.24 


-8.54 


-7.95 


-7.93 


-7.47 


-7.33 


DN Tau 




-8.89 


-7.50 


-8.46 


-8.72 


-8.39 


-8.89 


-7.50 


-8.26 


-7.56 


-6.93 


DO Tau 




-7.22 


-5.60 


-6.84 


-7.15 


-6.70 


-7.22 


-5.60 


-6.82 


-6.17 


-6.17 


DR Tau 






-5.10 




-6.66 


-5.88 


-6.66 


-5.10 


-8.16 


-6.17 


-6 


FT Tau 


















-8.56 


-7.29 


-6.79 


GG Tau 




-7.52 


-6.70 


-7.76 




-7.33 


-7.76 


-6.70 


-7.88 


-7.16 


-6.58 


GM Aur 




-8.13 


-7.60 


-8.02 


-8.11 


-7.97 


-8.13 


-7.60 


-8.5 


-7.37 


-7.37 


HL Tau 


















-6.53 


-o.yy 


c on 


IQ Tau 




-7.74 




-7.55 




-7.65 


-7.74 


-7.55 


-8.26 


-6.93 


-6.71 


IRAS 040164 


2610 
















-7.22 


-6.26 


-5.41 


IRAS 041694 


-2702 
















-7.76 


-6.73 


-6.12 


IRAS 042484 


2612 
















-10.41 


-7.57 


-6.89 


IRAS 043024 


-2247 
















-11.2 


-9.3 


-7.33 


IRAS 043254 


2402 
















-11.1 


-8.88 


-7.28 


IRAS 043614 


2547 
















-11.97 


-8.56 


-5.55 


IRAS 0436542535 
















-8.6 


-6.26 


-5.68 


IRAS 043684 


2557 
















-11.01 


-11.01 


-7.3 


L1551 IRS5 


















-9.29 


-5.34 


-4.93 


LkCa 15 








-8.87 




-8.87 


-8.87 


-8.87 


-8.55 


-8.33 


-7.16 


RY Tau 






-7.60 






-7.60 


-7.60 


-7.60 


-6.61 


-5.98 


-5.92 


UY Aur 




-8.31 


-6.60 


-7.18 




-7.36 


-8.31 


-6.60 


-8.94 


-8.57 


-6.59 



Note. — All disk accretion rates are in log [M Q / yr] . 

References. - V93: IValenti et all |l993l ). H95: lHartigan etal] (| 19951) . H98: iHartmann et al. 



(|l998h . M05: iMohantv et all (|2005l) 
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Table 5. Range of parameter values providing a good fit from SED fitting for additional 

disk parameters. 



Source Name 


^disk 


^>disk 


^j>disk 


y^disk 


odisk 
max 


odisk 
max 


/i(100AU) 


/i(100 AU) 




Min 


Max 


Min 


Max 


Min 


Max 


Min 


Max 




(-Rsub) 


(Rsub) 


(AU) 


(AU) 


(AU) 


(AU) 


(AU) 


(AU) 


AA Tau 


1.0 


6.4 


0.04 


0.51 


62.1 


168.9 


2.33 


6.67 


AB Aur 


1.0 


2.7 


0.48 


1.19 


83.1 


170.7 


2.90 


7.28 


BP Tau 


1.0 


6.4 


0.05 


0.51 


19.3 


585.8 


1.17 


9.26 


CI Tau 


1.0 


6.8 


0.08 


0.62 


77.8 


198.7 


1.17 


5.51 


CoKu Tau 1 


1.0 


1.0 


0.23 


0.23 


76.7 


76.7 


5.02 


5.02 


CY Tau 


1.0 


1.0 


0.04 


0.11 


97.5 


198.7 


1.17 


3.79 


DG Tau 


1.0 


6.1 


0.21 


2.50 


78.8 


2540.7 


2.69 


6.81 


DG Tau B 


1.0 


9.5 


0.07 


6.25 


2.3 


893.7 


2.76 


10.37 


DL Tau 


1.0 


1.0 


0.10 


0.11 


101.3 


198.7 


1.17 


3.16 


DM Tau 


130.2 


455.1 


3.37 


19.68 


50.4 


178.3 


3.53 


11.64 


DN Tau 


1.0 


6.4 


0.06 


0.51 


49.6 


163.1 


2.33 


5.77 


DO Tau 


1.0 


1.5 


0.16 


0.26 


103.7 


190.6 


2.42 


2.91 


DR Tau 


1.0 


5.8 


0.21 


1.39 


72.4 


158.5 


2.26 


7.61 


FT Tau 


1.0 


12.1 


0.04 


1.01 


24.8 


489.3 


1.17 


17.42 


GG Tau 


1.0 


6.8 


0.06 


0.62 


55.0 


133.0 


3.14 


12.47 


GJV1 Aur 


12.1 


265.8 


1.01 


oo no 

23.93 


69.3 


352.2 


2.30 


5.84 


HL Tau 


1.0 


2.9 


0.17 


0.43 


2.2 


39.0 


4.80 


19.30 


IQ Tau 


1.0 


6.4 


0.07 


0.51 


49.6 


415.8 


1.17 


5.77 


IRAS 04016+2610 


1.0 


2.9 


0.09 


2.46 


2.2 


3651.1 


3.35 


30.16 


IRAS 04169+2702 


1.0 


3.9 


0.07 


0.42 


3.6 


150.7 


5.28 


27.90 


IRAS 04248+2612 


1.0 


48.9 


0.03 


2.20 


27.0 


211.5 


2.74 


13.97 


IRAS 04302+2247 


77.9 


1700.5 


3.25 


58.85 


36.2 


1215.0 


3.31 


10.12 


IRAS 04325+2402 


95.6 


497.5 


4.29 


21.55 


22.8 


393.5 


3.69 


14.99 


IRAS 04361+2547 


1.0 


11.3 


0.07 


1.63 


2.2 


1296.0 


3.43 


24.36 


IRAS 04365+2535 


1.0 


9.3 


0.09 


1.60 


2.7 


176.7 


5.11 


15.54 


IRAS 04368+2557 


1.0 


2.5 


0.05 


0.24 


33.8 


1303.3 


3.18 


9.67 


L1551 IRS5 


1.0 


30.0 


0.17 


12.10 


1.2 


731.2 


2.81 


21.12 


LkCa 15 


1.0 


2.7 


0.04 


0.12 


62.1 


195.5 


4.42 


9.84 


RY Tau 


1.0 


4.7 


0.29 


1.38 


37.9 


124.1 


1.92 


5.76 


UY Aur 


1.0 


10.1 


0.27 


1.94 


51.5 


1159.7 


3.47 


10.44 
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Table 6. Range of parameter values providing a good fit from SED fitting for the stellar 

mass and total bolometric luminosity. 



Source Name M* M* Lboi ^boi 

Min Max Min Max 



(M ) (M ) (L ) (L ) 



AA Tau 


0. 


.20 


1 


.74 


0.44 


3.32 


AB Aur 


2. 


.57 


3 


.62 


42.64 


181.47 


BP Tau 


0. 


.22 


1 


.74 


0.59 


4.48 


CI Tau 


0. 


.35 


1 


.78 


1.60 


5.12 


CoKu Tau 1 


2. 


.68 


2 


.68 


11.44 


11.44 


CY Tau 


0, 


,57 


1 


.08 


0.40 


3.01 


DG Tau 


1, 


.10 


3 


.07 


9.63 


42.02 


DG Tau B 


0. 


,11 


3 


.70 


0.67 


249.65 


DL Tau 


0. 


,99 


1 


.29 


2.53 


3.01 


DM Tau 


0. 


,10 





.20 


0.16 


0.44 


DN Tau 


0. 


,33 


1 


.29 


0.82 


2.70 


DO Tau 

1 J V J ± Cl U 


n 

w. 


80 


I 


49 


5 40 


7 59 


DR Tau 


0. 


,63 


2 


.13 


7.53 


14.36 


FT Tau 


0, 


,11 


2 


.03 


0.36 


7.01 


GG Tau 


0. 


,28 


1 


.63 


0.82 


3.01 


GM Aur 


0. 


,28 


1 


.62 


0.65 


3.16 


HL Tau 


0. 


,10 





.91 


3.46 


15.05 


IQ Tau 


0. 


,33 


1 


.94 


1.04 


5.16 


IRAS 04016+2610 


0. 


.10 


3 


.73 


2.00 


249.21 


IRAS 04169+2702 


0. 


,11 





.90 


0.73 


5.87 


IRAS 04248+2612 


0. 


,10 


2 


.84 


0.27 


11.30 


IRAS 04302+2247 


0. 


,10 





.25 


0.15 


0.67 


IRAS 04325+2402 


0. 


.10 





.41 


0.46 


1.54 


IRAS 04361+2547 


0, 


,11 





.85 


1.20 


9.37 


IRAS 04365+2535 


0. 


,11 


1 


.25 


2.00 


14.43 


IRAS 04368+2557 


0, 


.15 





.59 


0.64 


3.84 


L1551 IRS5 


0, 


,21 


4 


.70 


6.88 


72.40 


LkCa 15 


0. 


,18 





.91 


0.42 


1.39 


RY Tau 


1. 


,64 


3 


.34 


18.75 


40.65 


UY Aur 


0, 


,37 


2 


.04 


2.11 


14.22 
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Table 7. The main parameters for the best-fitting models to IRAS 04368+2557 



Parameter 


IRAC data only 


Full SED 


Stellar Mass 


4.07 M Q 


1.46 M 


Stellar Radius 


21.57 R Q 


8.21 R Q 


Stellar Temperature 


4360 K 


4260 K 


Envelope accretion rate 


2.63xlO" 4 M /yr 


1.37xl0" 4 M o /yr 


Envelope outer radius 


9120 AU 


16200 AU 


Cavity opening angle 


16° 


43° 


Viewing angle 


75° 


81° 


Bolometric Luminosity 


155 L & 


20 L Q 



Note. - These are the main parameters of the best-fit model 
SEDs to the SED of IRAS 04368+2557, first fitting only the IRAC 
data simulaneously in six apertures, then fitting the full SED in- 
cluding the multi-aperture IRAC data. 
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A. The convolution of model SEDs with broadband filters 

This appendix describes the exact procedure used to obtain monochromatic fluxes 
through broadband filters for our models. 

In the following, we define the true spectrum of a source or model SED to be F v [actual]. 
In general, when a broadband flux is measured through a filter, we have no knowledge of 
the true underlying spectrum, only the integrated flux over the filter. Therefore, to quote a 
monochromatic flux F UQ [quoted] at a frequency u , one usually makes an assumption about 
the spectrum of the source. We call this spectrum F„ [assumed] . Note the difference between 
F v [assumed] and F VQ [quoted] : F VQ [quoted] is the value of F„ [assumed] at z/ . 

We make the same assumption for our model SEDs as is made for the data taken in the 
different filters. For example, fluxes from the IRAC pipeline are quoted using ^[assumed] oc 
\jv. That is to say, what would the flux at u be, if the spectrum of the source was 
proportional to 1/u, such that the observed integrated flux was identical to what is actually 
observed? Since we are making the same assumptions as used for the observed fluxes, a 
direct comparison between the quoted model fluxes and the quoted observed fluxes can be 
made without the need for any color-correction. 



A.l. Spitzer - IRAC 



The IRAC mo nochromatic fluxes assume ^[assumed] oc or vF v [assumed] oc const 
faeach et al.ll2005h . 



The total electron "count" detected through the filter with response R{v) (in e /photon) 

is 

E = [ F4 f Ual] R{v) du= f F - [aS f Umed] R(u) dv. (Al) 
J hv J hv 

Now vF v [assumed] oc const = vqF Vo [quoted] , so 

v F UQ [quoted] f R{y) dv = [ F ^ tuaI ] R (y) dv ( A2 ) 
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which after rearranging gives 

= / factual] (, H flfr) A/ 
° L J J (v /vf R(v) dv 

The values of i_(v) are taken from the Spitzer Science Center websita_i The values of 
Uq are given by v = c /A n where An = 3 .550, 4.493, 5.731, and 7.872/im are the nominal 



wavelengths for IRAC (IReach et al.ll2005l ) 



A.2. Spitzer - MIPS 

The MIPS monochromatic fluxes assume a T=10, 000K blackbody spectrum, i.e very 
close to F„ [assumed] oc v 2 , or F u [assumed] jv 2 oc const (MIPS data handbook^). 

The total electron "count" through the filter with response R v {v) = R(v) / v (in e~/unit energy) 

is 

E = J F u [actual] R v (v) dv = J F„ [assumed] R v {v) dv. (A4) 
Now F,, [assumed] /z/ 2 oc const = F VQ [quoted] /vq, so 

F UQ [quoted] /v 2 J v 2 R v {y) dv = J F„[actual] R v {v) dv (A5) 
which after rearranging gives 

f F v [actual] R v (v) dv , , , 

F uo [quoted] = J —f- 1 v \> (A6) 
J [v/vo) R v (v) dv 

The values of R v {y) are taken from the Spitzer Science Center website,]- The values of v 
are given by vq = c/An where Ao = 23.68, 71.42, and 155.9/xm are the effective wavelengths 
for MIPS. 



5 http:/ /ssc. spitzer. caltech.edu/irac/spectral_response. html 

6 http:/ /ssc. spitzer. caltech.edu/mips/dh/ 

7 http://ssc. spitzer. caltech.edu/mips/spectral_response. html 
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A.3. IRAS 

The IRAS monochromatic fluxes use ^[assumed] oc or vF v [assumed] oc const 
(IRAS Explanatory Supplement - Section VI.C.3^|). In the same way as for IRAC, we have 

F U0 [quoted] = H^^M g&j ± (A7) 

However, the relative system response listed in the IRAS documentation is in electrons 
per unit energy (as for MIPS). Therefore, writing R v (v) = R{y)/v, we get 

„ r f -FJactuaf] R v {v) dv ... 

F U0 quoted = j u J ' A8 

J (z/ /f) ^M^) «^ 

The values of R u {v) are taken from the IRAS documentation (The 'Relative System 
Response' in the IRAS Explanatory Supplement - Table II.C.Eq). The values of u are given 
by pq = cj A where A = 12, 25, 60, and lOOfim. 



A.4. UBVRI photometry 



The UBVRI observations from lKenyon fc Hartmannl (119951 ) are originally from lHerbst et al 



(119941 ). and were made in the J ohnson/ Cousi ns UBVRI system. The transmission curves 
for these bands were taken from iBesselll (jl990j). Since flux densities are not commonly used 
at optical wavelengths, it was not clear what spectrum to assume in order to derive the 
monochromatic fluxes (although the calibration is usually done using the spectrum of Vega). 
However, we have found that the differences arising from various assumptions do not change 
the resulting fluxes by more than a few %, which is much smaller than the ±25% uncertain- 
ties we imposed on the observed UBVRI fluxes. Therefore, the choice of the assumption is 
unimportant for this work. We choose to assume a flat spectrum, i.e. ^[assumed] oc 1/V, 
or vF v [assumed] cx const. Therefore, as for IRAC, we have: 



F Uo [quoted] 



f F v [actual] {y^jv) R(y) dv 



(A9) 



8 http:/ /irsa. ipac.caltech.edu/IRASdocs/exp.sup/ch6/C3. html 
9 http://irsa. ipac.caltech.edu/IRASdocs/exp.sup/ch2/tabC5. html 
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The values of u are given by u = c/X where Aq = 0.36, 0.44, 0.55, 0.64, and 0.79/wn. 



A. 5. JHK/2MASS 

For all the JHK flu xes we computed the m onochromatic fluxes using the method for the 
2MASS all-sky survey Jskrutskie et al.lbood^ 

The 2MASS is ophotal fluxes are computed 
using the relative system response R v from ICohen et al.l (120031 ) . and are given by: 



^[quoted] 



f F u [actual] R v {y) di 
Az/ isn 



(A10) 



The isophotal bandwidths Az/ iso for the three bands are listed in I Cohen et al.l (120031 ). 
The isophotal wavelengths are Ag° = 1.235, 1.662, and 2.159/im. 



A.6. SHARC II, SCUBA, and CSO observations 

The SHARC II and SCUBA instruments are calibrated on planets, whose radiation fol- 
lows the Rayleigh- Jeans tail of a blackbody curve at sub-mm wavelengths. Therefore, we 
take ^[assumed] oc z/ 2 , or F v [assumed] jv 2 oc const. 

As for MIPS, we have 



F„ [quoted] 



f F u [actual] R v {y) dv 



(All) 



/ {v/vof Ru(v) dv 

The values of R u {v) are taken from lDowell et al.l J2003h for the SHARC II observations, 
and the SCUBA websitfi for the SCUBA 450WB and 850WB observations. For the CSO 
observations we used gaussians centered at 624 and 769 /im with FWHM 67 and 190 ^m re- 
spectively. We convolved the S HARC II and CSO filters with the atmospheric transmission 
curve used in lDowell et al.l (120031 ). and the SCUBA filters with the atmospheric transmission 
curve given with the filter profiles on the SCUBA website. 



3 http://www .jach.hawaii.edu/JCMT/continuum/background/background. html 



-47- 



The values of u are given by v = c/A where A = 350, 443, and 863/xm for SHARC II 
350 fim and SCUBA 450WB and 850WB respectively, and A = 624 and 769^m for the CSO 
observations. 
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B. The data for the 30 Taurus- Auriga sources 
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Table 8. Optical and Near-IR data for the 30 Taurus- Auriga sources. 



Source Name 


U 


B 




V 


R 


I 




J 




II 


K 


References 


A A Tau 


3.60 


8.97 


26 


24 


48. 11 


104.42 


271 


95 


389 


07 


404.70 


1 . 5 




± 1.61 


± 3.62 


± 10 


20 


± 18.80 


± 40.92 


± 7 


76 


± 13 


28 


rh 14.53 




AB Aur 


2361.72 


5547.50 


5723 


75 


4837.74 


4797.35 


6731 


14 


9671 


64 


13549.72 


1 . 5 




± 106.45 


+_ 144.47 


± 105 


42 


± 125.98 


± 152.98 


+_ 111 


57 


± 178 


13 


+_ 199.64 




BP Tau 


9.87 


20.89 


51 


93 


88.44 


175.02 


306 


85 


546 


39 


547.94 


1 . 5 




± 2.01 


± 2.83 


± 4 


76 


± 8.73 


± 17.93 


± 13 


18 


± 15 


60 


± 15.64 




CI Tau 


3.57 


7.05 


18 


95 


38.77 


95.68 


260 


91 


438 


43 


504.35 


1 . 5 




± 1.29 


± 2.07 


± 3 


60 


± 7.36 


± 18.17 


± 6 


73 


± 17 


36 


rh 13.93 




CoKu Tau 1 










1.11 


3.41 


11 


38 


25 


98 


27.18 


3 . 5 












± 0.00 


rh 0.00 


± 


32 


± 


75 


± 0.61 




CY Tau 


3.01 


5.38 


17 


20 


33.92 


89.30 


197 


01 


268 


35 


249.07 


1 . 5 




± 0.64 


± 0.96 


± 1 


58 


± 3.17 


± 8.70 


± 5 


81 


± 9 


64 


± 7.34 




DG Tau 


10.93 


18.04 


42 


55 


77.27 


169.33 


532 


22 


834 


64 


1064.46 


1,5 




± 3.60 


± 5.79 


± 13 


37 


± 24.42 


± 53.97 


± 53 


22 


± 83 


46 


+_ 106.45 




DG Tau B 














1 

± 


95 

19 


5 

± 


91 
59 


14.24 
± 1.42 


5 


DL Tau 


5.14 


8.41 


21 


80 


55.52 


120.44 


229 


13 


352 


13 


438.06 


1 . 5 




± 1.29 


± 2.00 


± 4 


33 


± 14.65 


+_ 31.79 


± 7 


39 


± 10 


70 


± 12.91 




UJVL lau 


2.57 


3.99 


9 


54 


20.30 


52.47 


100 


09 


128 


20 


103.83 


1.5 




± 0.69 


± 0.80 


± 1 


56 


± 3.62 


± 9.38 


± 3 


03 


± 3 


42 


± 3.16 




DN Tau 


2.96 


12.29 


41 


19 


74.77 


168.34 


345 


53 


474 


13 


415.27 


1 . 5 




± 0.66 


± 1.06 


± 3 


03 


± 5.66 


± 13.56 


± 17 


82 


± 21 


83 


± 12.24 




DO Tau 


1.91 


3.84 


12 


75 


35.84 


101.15 


261 


39 


520 


84 


801.55 


1.5 




± 1.01 


± 1.68 


± 4 


87 


± 15.76 


± 45.33 


± 2 


65 


± 18 


23 


± 25.10 




DR Tau 


44.27 


56.20 


98 


50 


140.72 


254.35 


470 


97 


761 


20 


1192. 15 


1,5 




± 21.70 


± 26.95 


± 46 


64 


± 66.66 


± 121.11 


± 16 


69 


± 41 


36 


± 40.62 




FT Tau 




2.32 


5 


03 


22.56 


75.30 


143 


12 


239 


39 


255.35 


2.5 






± 0.46 


± 


00 


± 4.48 


± 0.00 


± 4 


22 


± 7 


27 


± 7.29 




GG Tau 


6.55 


14.88 


49 


44 


94.85 


193.02 


547 


13 


757 


71 


764.07 


1.5 




± 0.80 


± 0.80 


± 2 


27 


± 4.70 


± 11.90 


± 27 


71 


± 25 


12 


± 30.25 




GM Aur 


6.93 


21.28 


58 


29 


95.61 




299 


56 


379 


40 


322.06 


1.5 




± 0.00 


± 0.00 


± 


00 


± 0.00 




± 9 


93 


± 9 


08 


± 7.12 




HL Tau 


0.49 


1 .49 


4 


95 


10.12 


25.21 


89 


72 


219 


74 


724.32 


1.5 




± 0.11 


± 0.23 


± 


72 


± 1.49 


± 3.71 


± 3 


47 


± 9 


31 


± 11.34 




IQ Tau 


1.59 


5.20 


1 7 


00 


39.34 


103.57 


270 


70 


454 


46 


522.79 


1.5 




± 0.94 


± 2.93 


± 6 


24 


± 15.46 


± 40.78 


± 9 


72 


± 14 


23 


± 19.26 




IKAo U4Ulu + z610 









03 


0.06 


0.18 


5 


98 


22 


24 


123.62 


1.4 








± 


00 


± 0.00 


± 0.00 


± 1 


84 


± 3 


30 


± 3.41 




IRAS 04169+2702 
















± 


34 
17 


3 

± 


01 
10 


21.68 
± 0.40 


4 


IRAS 04248+2612 














12 

± 2 


73 
12 


28 

± 4 


83 
58 


37.09 
± 5.74 


4 


IRAS 04302+2247 














3 

± 


33 
00 


11 
± 


76 
00 


24.88 
± 0.00 


4 


IRAS 04325+2402 


















6 

± 


34 
00 


25.82 
± 0.00 


4 


IRAS 04361+2547 
















± 


61 
13 


7 

± 


28 
96 


40.28 
± 2.96 


4 


IRAS 04365+2535 
















± 


20 
05 


4 

± 1 


64 
10 


38.51 
± 7.98 


4 


IRAS 04368+2557 














5 

± 2 


05 
52 


31 

± 15 


09 

55 


75.39 
± 37.69 


5 


L1551 IRS5 














10 

± 5 


88 
48 


47 

± 23 


06 
21 


133.53 
± 27.67 


4 


LkCa 15 


2.75 


13.59 


43 


87 


73.39 


136.38 


270 


45 


375 


58 


357.05 


1,5 




± 0.97 


± 4.73 


± 15 


13 


± 25.57 


± 48.12 


± 5 


73 


± 10 


03 


± 10.52 




RY Tau 


42.79 


146.36 


345 


54 


507.61 


891.41 


2255 


74 


3797 


52 


4746.03 


1,5 




± 18.92 


± 64.16 


± 150 


72 


± 221.69 


± 390.17 


± 66 


47 


± 220 


31 


± 122.37 




UY Aur 


3.88 


9.26 


27 


51 


57.55 


141.73 


354 


89 


652 


68 


866.82 


1.5 




± 1.83 


± 4.35 


± 12 


86 


± 27.20 


± 68.36 


± 10 


46 


± 13 


82 


± 25.54 





References. — 1: UBVRI data from Kcnyon & Hartmann 1 1 995^ . 2: B VRI data from the USNO B -1.0 catalog, 3: average BVRI data 
from the USNO B-1.0 catalog and the NOMAD catalog, 4: JHK data from lKenvon &: Hartmannl il995f) . and 5: JHK data from the 2MASS 
all-sky survey (the fluxes for DG Tau, DG Tau B, and IRAS 04368 + 2557 were measured using aperture photometry) 



Note. — Fluxes are in mjy. Values in italics indicate that these are not used when fitting the observed SEDs as higher quality values 
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arc available 
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Table 9. Mid-IR data for the 30 Taurus-Auriga sources. 



Source Name 


L 


M 


N 


Q 


IRAC 3.6 (im 


IRAC 4.5 


IRAC 5.8 /im 


IRAC 8.0 /im 


References 


J\J\ X au 


126 29 


297 0# 


456 14 




.'-i 70 


5 5 


352 


24 


331 


89 


3 5 5 


82 


I 2 




i 77 22 




-j- go 






(>4 


_j_ 12 


97 


zL 12 


v> 2 


i 9 


83 




AB Aur 


1 97^.8 ^7 
± z / JO . o / 

zn ^uzz.o^ 


1 1 SH9 fi9 
-1- ^007 91 


90001 p>i 

4- / R9K 9R 

ZTZ /(OZJ.iiO 


O 1 ±0 . IH 

A- 1 QR&K P I 
zn ijflij.i^ 


















1 


BP Tau 


5S£> 75 




500 S6 


6'JS 7J 






281 


25 


2 3 1 


23 


3 3 5 


80 


1 4 








+ 109.58 








i 3 








i 16 






CI Tau 


5i5 Oi 




673 23 




()4() 


10 


5 5 8 


67 


487 


90 


572 


97 


1 4 








_|_ ^ tjg 




-j- g 


08 


i 34 


1 1 


zL 11 


14 


i 25 


3 5 




CoKu Tau 1 


23.96 








2 1 


72 


49 


5 5 


102 


98 


3 1 5 


60 


1.2 




i j yg 








_J_ Q 


s Ji 


i 2 


" S 


zb 1 




zb 5 






CY Tau 


Si.? ^4" 






















9r 


I 2 












_J_ g 


C2 


zb 9 


4() 


zb 4 


11 


zb 2 


49 




DC lau 


£.0 U d. O U 


1800 51 


/ {){) :JU 






HO 


2337 


04 


2097 


00 


3307 


00 






i 567 06 


-A- 1 1 19 77 
zrz -/ .£ cs. /4 


4- 1 fiQfi 01 
zrz J eye. t/J 


zn .)u4^-4 ' 


dz 56 


63 


zL 153 


35 


+ 254 


70 


± 111 


03 




TIC 1 Tan R 


98 41 


00 


00 


00 


( >0 


(>2 


287 


88 


571 


47 


830 


98 


1 2 




i 00 


i 00 


i 00 


-J- Q QQ 


=b 5 


00 


i 15 


88 


± 21 


04 


zrz 22 


93 




FIT Tan 


60 4 38 








4S4 


00 


533 


14 


5 2 5 


85 


590 


51 


1 2 




i 70 


+ 65 37 


+_ iS7 81 




4. 27 


82 




(>3 


i 14 


52 


± 16 


31 




DM Tau 


1 05.50 




Ann Qfi 




55 


72 


37 


80 


23 


24 


22 


13 


1 4 




+ 86.60 




4- 1 90 09 




± 


37 


+_ 1 


49 


± 3 


52 


± 2 


70 




DN Tau 


812.43 


324 . 93 


261 . 63 




3{). r i 


50 


202 


55 


242 


00 


2 5 ;-; 


06 


1.2 




+ 72.93 


+_ 41.55 


+ 42.33 




± 14 


05 


-j— 14 


49 


± 8 


93 


± 


99 




DO Tau 


951 62 


19 lf> 8/ 


1481.31 


9/80 J 1 


965 


42 


988 


57 


95 


89 


1026 


19 


1.2 




4- 1 RO RR 
zn luuiOO 


zn ± &u . CJ4 


-A- 91R '1R 
zjz <ouo .00 


i 204.91 


4. 26 


(>(> 


+ 45 


48 


-]_ 26 


43 


zb 28 


34 




FIR Tun 


231 5 .37 


2839 78 


1O0R QR 




1 8 . r ) 8 


72 


189(3 


13 


'■'004 


98 


1736 


27 


1 4 




3Z OUU. Z>4 


-A- 1 1 10 QQ 
zn 1 iou . c?f 


+ 1115 34 




4. 204 


37 


i 151 


04 


i 304 


25 


i 216 


95 




FT Tau 


SOS ^0 








24') 


5 3 


'• > 48 


80 


205 


02 


271 


85 


1 4 




i 00 








4. 2 








zb 1 




i 10 






Kj ± a li 




722 98 


945 58 


/ 90 / 9 


671 


(ifi 


5 7 


00 


472 


32 


559 


22 


1 4 




i 50 07 




zfc 248.12 


± 46.18 


4. 33 


22 


i 10 


93 


_j_ ig 


54 


-j- 24 


92 




1 VI rt. U I 


i 61 34 




461 98 




171 


f j 3 


120 


87 


95 


66 


102 


42 


1 3 




i 25.08 




+ 67.34 




± 2 


2 1 


zrz 


70 


zb 


7!) 


— |— 2 


36 




TZTT Tan 


1 R1 1 OR 
1 1 4 ■ l/o 


I OfiH on 
4 (7 u 1 .:JU 


fi 1 %n or 


9 1 R1 / f\9 

eC lo 14 .O £ 


3142 


80 


4285 


42 


5119 


98 


4470 


30 


1 4 




_i_ ijflo an 


■A- 1 101 81 

ZTZ -Z JtZ-i . O-i 


_i_ 9'yyo 1 1 
zrz iC 1 1 -44 


-A- 101 r .R ,y >8 
zn 1 (>/j i>o . -0 






i 41 














TO Tan 


^iS 4" 6 




^05 47 




48r> 


54 


450 


15 


384 


07 


377 


29 


1 2 






' ' ' 






it 8 




+. 3 




zb 8 




zb 2 






tdaq nzim fij_9fii n 

lr\.^Vo IMU1 1> f ZU J. U 
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1 4 5 3 


51 


1708 
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31 
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^ 227 IS 




+ 9S><? 8 
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195 


72 


300 


10 


333 


10 


417 


53 


1 4 




± 2 89 




± i0^0 




± 9 


27 


± 9 


99 


± 3 


49 


± 14 


07 




IRAS 04248+2612 


28.44 


' . . 


256.55 


96.9.00 


65 


22 


67 


36 


78 


44 


109 


88 


1,2 




± 8.99 




± 39.68 


± 149.88 


± 5 


98 


± 6 


18 


zb 7 


19 


zb 10 


08 




IRAS 04302+2247 


21.40 








46 


25 


39 


35 


27 


43 


12 


68 


1,4 




± 0.00 








± 


76 


± 


46 


zb 


03 


zb 1 


17 




IRAS 04325+2402 










74 


88 


77 


34 


58 


96 


38 


81 


2 












± 6 


87 


± 7 


0!) 


± 5 


41 


± 3 


50 




IRAS 04361+2547 


1 08. 71 




1059.40 




264 


48 


353 


50 


415 


47 


880 


85 


1,4 




± 28.49 


zb J 7. SO 


± 77. S5 




± 24 


20 


± 32 


42 


zb 38 


10 


zb 80 


7!) 




IRAS 04365+2535 


139.33 


S76.5£ 


650.91 


5280.44 


374 


30 


892 


98 


1421 


88 


1469 


69 


1,2 




± 28.87 


zb 57.50 


± 


± 1094.05 


± 20 


05 


± 32 


88 


± 39 


27 


zb 40 


59 




IRAS 04368+2557 










166 


87 


221 


00 


264 


57 


82 


58 


2 












± 15 


30 


± 20 


27 


zb 24 


27 


zb 7 


57 




L1551 IRS5 


447.41 


J^SS.O^ 


5266.^7 


25042.07 


598 


89 


1375 


48 


3105 


00 


4908 


54 


1.4 




± 179.84 


zb 814-78 


± 1091.15 


± 5188. 44 


± 5 


83 


± 2 


70 


± 34 


90 


zb 41 


39 




LkCa 15 


253.14 




201.46 




267 


79 


207 


04 


142 


68 


163 


2!) 


1.4 




± 16.77 




± 7.^2 




± 1 


53 


± 2 


62 


± 3 


77 


zb 1 


00 




RY Tau 


5677.26 


7554.89 


15922.10 




5749 


96 


5298 


11 


4403 


50 


5829 


60 


1.4 




± 1057.69 


=b 2110.88 


± 5629.73 




± 72 


53 


± 57 


59 


± 157 


21 


± 161 


56 




UY Aur 


1398.63 




2809.14 


4939.96 


1072 


04 






1106 


39 


1713 


42 


1,3 




± 258.84 




± ^09.52 


± 720.1 6 


± 17 


77 






zb 8 


15 


zb 4 


73 





References. — 1: LMNQ data from Kcnyon & Hartmann (1995), 2: IRAC data from Hartmann ct al. (2005), 3: IRAC data from Luhman ct al. (2006), 4: IRAC 
data measured using data retrieved from the Spitzcr Space Telescope Archive 

Note. — Fluxes arc in mjy. Values in italics indicate that these arc not used when fitting the observed SEDs as higher quality values are available 
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Table 10. Far-IR dat 


a for the 30 Taurus- Aurig 


- a sources. 




Source Name 


MIPS 24 


MIPS 70 


IRAS 12 


IRAS 25 


IRAS 60 


IRAS 100 


References 


AA Tau 


502.83 


950.05 


430.00 


610.00 


1230.00 


QQfin nn 
a&ifU- UU 


1.2 




± 2.80 


± 95.00 


± 33.00 


± 32.00 


± 96.00 


zt 601 . 00 




AB Aur 






28950.00 


49780.00 


115680.00 


1 1 A A *~7(\ AA 

1 144 /U.UU 


2 








± 63.00 


± 107.00 


± 72.00 


1 1 rtrt nn 

± lyob.uo 




BP Tau 






450.00 


590.00 


440.00 


non aa 

y^u.uu 


2 








± 35.00 


± 36.00 


± 2.00 


zt 226.00 




CI Tau 


956.95 




780.00 


1300.00 


2150.00 


< zo4U.U0 


1,2,3 




± 5.70 




± 25.00 


± 46.00 


± 71.00 






CoKu Tau 1 


3453.61 




1180.00 


2740.00 


< 7970.00 


^ "71 nnn nn 

<c * luyu.uu 


1,2,3 




± 38.10 




± 26.00 


± 63.00 








CY Tau 


130.48 


98.86 


190.00 


270.00 


140.00 




1,2 




± 1.38 


± 9.89 


± 45.00 


± 35.00 


± 90.00 






DG Tau 


9332.72 


9824.49 










1 




± 854.50 


± 982.45 












DG Tau B 


4641.72 


7836.02 










1 




± 99.20 


± 783.60 












DL Tau 


906.17 


888.20 


970.00 


1320.00 


1390.00 


<?>o on nn 
3S&3U. 00 


1,2 




± 4.80 


± 88.82 


± 34.00 


± 45.00 


± 84.00 


i n 1 n nn 

± 94 -J- 00 




DM Tau 






< 270.00 


350.00 


830.00 


^- 701 A A A 

< (zlU.UU 


2,3 










± 33.00 


± 104.00 






DN Tau 


410.03 


423.01 


350.00 


600.00 


650.00 


< 5787.00 


1,2 




± 2.60 


± 42.30 


± 30.00 


± 53.00 


± 93.00 






DO Tau 


3115.98 


2745.71 


1880.00 


4070.00 


6330.00 


c c n nn 


1,2 




± 32.30 


± 274.57 


± 34.00 


± 41.00 


± 106.00 


1 1 too nn 
zt I400.OO 




DR Tau 






3160.00 


4300.00 


5510.00 


cfiufi nn 


2 








± 29.00 


± 51.00 


± 44.00 


1 11 AA 

zt 1144.UU 




FT Tau 






360.00 


570.00 


820.00 


1 onn nn 
loUU.UU 










± 45.00 


± 37.00 


± 41.00 


I 1 Q "7 A A 




GG Tau 






1200.00 


1720.00 


3120.00 


r,ion nn 
O4OU.U0 


2 








± 33.00 


± 39.00 


± 71.00 


1 OCC A A 




GM Aur 


747.62 


1877.76 


250.00 


1070.00 


3080.00 


i i n nn 

0440 -00 


1,2 




± 6.55 


± 17.80 


± 31.00 


± 42.00 


± 112.00 


1 1 nan nn 




HL Tau 






9740.00 


31180.00 


76260.00 


'j'mrzn nn 
/ f yoU.UU 










± 31.00 


± 73.00 


± 810.00 


I 1 A A A A 

zh _L.iy4.UU 




IQ Tau 


500.09 


579.59 


500.00 


650.00 


770.00 


^ 1 onn nn 


1,2,3 




± 3.00 


± 57.96 


± 29.00 


± 30.00 


± 43.00 






IRAS 04016+2610 


11992.36 




3640.00 


15810.00 


48790.00 


ccKfin nn 

oooyu.uu 


1,2 




± 257.11 




± 109.20 


± 632.40 


± 3903.20 


I Tn 




IRAS 04169+2702 


4519.14 


11745.45 


750.00 


5210.00 


17000.00 


f/v/ f n nn 

1 14.00. 00 


1,2 




± 821.84 


± 1174.55 


± 52.50 


± 312.60 


± 1 700. 00 


1 ^> fj) /? r J on 
zt HZv>). <i0 




IRAS 04248+2612 


835.33 


3640.03 


< 360.00 


1330.00 


4620.00 




1.2 




± 4.80 


± 364.00 




± 106.40 


± 415. 80 


1 000 in 
zt 000.40 




IRAS 04302+2247 


241.09 


4774.75 


< 250.00 


440.00 


6400.00 


9430. 00 


1.2 




± 1.80 


± 477.48 




± 83.60 


± 6^0.00 


zt 1131.60 




IRAS 04325+2402 


1860.53 


8042.05 


< 250.00 


2100.00 


12860.00 


22350.00 


1.2 




± 11.60 


± 804.21 




± 168.00 


± 1157.40 


± 3576.00 




IRAS 04361+2547 






1820.00 


18870.00 


44750.00 


35430.00 


2 








± 109.20 


± 1132.20 


± 5370.00 


± 4251.60 




IRAS 04365+2535 


6118.02 




1190.00 


8620.00 


36010.00 


39250.00 


1,2 




± 399.90 




± 142.80 


± 517.20 


± 4321.20 


± 5887.50 




IRAS 04368+2557 


542.92 




< 250.00 


740.00 


17770.00 


73260.00 


1,2 




± 3.20 






± 66.60 


± 1599.30 


zt 11721.60 




L1551 IRS5 






10040.00 


106200.00 


372900.00 


457900.00 


2 








± 502.00 


± 4248.00 


± 18645.00 


zt 59527.00 




LkCa 15 






270.00 


390.00 


1500.00 


1610.00 


2 








± 30.00 


± 33.00 


± 48.00 


zt 200.00 




RY Tau 


17862.32 


9632.66 


17740.00 


26480.00 


18910.00 


13500.00 


1.2 




± 4424.47 


± 963.27 


± 27.00 


± 54.00 


± 66.00 


zt 2501.00 




UY Aur 






3710.00 


6870.00 


7580.00 


9400.00 


2 








± 43.00 


± 39.00 


± 87.00 


± 698.00 





Re ferences. — 1: MIPS d ata measured using data retrieved from the Spitzcr Space Telescope Archive, 2: IRAS data 
from IWeaver &: Jone 1 <1992tl . 3 : IRAS data from IRAS point-source catalog 



Note. — Fluxes arc in mjy. Values in italics indicate that these arc not used when fitting the observed SEDs as higher 
quality values arc available 
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Table 11. Sub-mm data for the 30 Taurus- Auriga sources. 



Source Name 


SHARC 350 iim 


SCUBA 450 pm 


SCUBA 850 /J,m 


CSO 624 fim 


CSO 729 (im 


References 


AA Tau 


825.00 


415.00 


144.00 




310.00 


1.2 




± 50.00 


± 84.00 


± 5.00 




± 60.00 




AB Aur 


8930.00 


3820.00 


359.00 






1 




± 1410.00 


± 570.00 


± 67.00 








BP Tau 




< 456.00 


130.00 






1 








± 7.00 








CI Tau 


1725.00 


846.00 


324.00 


1300.00 


850.00 


1.2 




± 55.00 


± 89.00 


± 6.00 


± 210.00 


± 150.00 




CoKu Tau 1 




< 522.00 


35.00 






1 








± 7.00 








CY Tau 


< 1839.00 


< 210.00 


140.00 




240.00 


1.2 








± 5.00 




± 40.00 




DG Tau 


< 5173.00 


< 3950.00 


< 1100.00 


< 1210.00 


< 860.00 


1.2 


DG Tau B 


< 5173.00 


< 3950.00 


< 1100.00 


< 1210.00 


< 860.00 


1.2 


DL Tau 


1390.00 


1280.00 


440.00 


880.00 


530.00 


1.2 




± 180.00 


± 170.00 


± 40.00 


± 140.00 


± 90.00 




DM Tau 


1077.00 




237.00 


390.00 


470.00 


1.2 




± 49.00 




± 12.00 


± 130.00 


± 80.00 




DN Tau 


615.00 


< 703.00 


201.00 




380.00 


1.2 




± 64.00 




± 7.00 




± 80.00 




DO Tau 




734.00 


258.00 


700.00 


510.00 


1.2 






± 50.00 


± 42.00 


± 100.00 


± 100.00 




DR Tau 




2380.00 


533.00 




400.00 


1.2 






± 172.00 


± 7.00 




± 80.00 




FT Tau 


1106.00 


437.00 


121.00 


260.00 


250.00 


1.2 




± 82.00 


± 56.00 


± 5.00 


± 100.00 


± 50.00 




GG Tau 


6528.00 


2726.00 


1255.00 


1370.00 


1250.00 


1.2 




± 153.00 


± 250.00 


± 57.00 


± 170.00 


± 8.00 




GM Aur 


3419.00 






1340.00 


850.00 


1,2 




± 133.00 






± 330.00 


± 90.00 




HL Tau 


23888.00 


10400.00 


2360.00 


5450.00 


3200.00 


1,2 




± 149.00 


± 1400.00 


± 90.00 


± 290.00 


± 100.00 




IQ Tau 




425.00 


178.00 






1 






± 26.00 


± 3.00 








IRAS 04016+2610 


12477.00 










1 




± 193.00 












IRAS 04169+2702 


7344.00 










1 




± 152.00 












IRAS 04248+2612 


1178.00 










1 




± 30.00 












IRAS 04302+2247 


2869.00 










1 




± 21.00 












IRAS 04325+2402 




606.00 


186.00 






1 






± 185.00 


± 11.00 








IRAS 04361+2547 




1302.00 


275.00 






1 






± 168.00 


± 8.00 








IRAS 04365+2535 




2928.00 


622.00 






1 






± 230.00 


± 13.00 








IRAS 04368+2557 




2849.00 


895.00 






1 






± 222.00 


± 11.00 








L1551 IRS5 


100423.00 










1 




± 812.00 












LkCa 15 


1235.00 




428.00 






1 




± 80.00 




± 11.00 








RY Tau 


2439.00 


1920.00 


560.00 






1 




± 330.00 


± 160.00 


± 30.00 








UY Aur 


542.00 


< 523.00 


102.00 






1 




± 77.00 




± 6.00 









Re ferences. — 1: SHARC II and SCUBA data from the compilation presented in lAndrews &: Williams] i2005T) . 2: CSO data 
from lBec kwith & Sargent <199lft 

Note. — Fluxes are in mjy. Values in italics indicate that these arc not used when fitting the observed SEDs as higher quality 
values arc available 
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Table 12. Apertures assumed for the SED fitting. 



a 

o o 1 1 rc c ~N amc 


UBVRI 




IRAC 


IVTTP^ OA urn 


Mil O 1 U ^1111 


TRAC io o. nc ,, m 


TR»c an e inn ,, m 


su b- nirn 


AA Tau 


5 


3 


5 


1(1 


120 


00 




30 


AB Aur 


5 


3 15 








00 


120 


30 


BP Tau 


5 


3 


5 






00 


120 


30 


CI Tau 


5 


3 


5 


10 




(iO 


120 


30 


CoKu Tau 1 


5 


3 


5 


10 




(iO 


120 


30 


CY Tau 


5 


3 


5 


10 


20 


(iO 




30 


DG Tau 


5 


10 


5 


10 


20 






30 


DG Tau B 




10 


5 


10 


20 






30 


DL Tau 


5 


3 


5 


10 


20 


(iO 




MO 


DM Tau 


5 


3 


5 






00 


120 


30 


DN Tau 


5 


3 


5 


10 


20 


(iO 




30 


DO Tau 


5 


3 


5 


10 


20 


(iO 




30 


DR Tau 


5 


3 


5 






00 


120 


30 


FT Tau 


5 


3 


5 






00 


120 


30 


GG Tau 


5 


3 


5 






00 


120 


30 


GM Aur 


5 


3 


5 


10 


20 


(iO 




30 


HL Tau 


5 


3 


5 






(iO 


120 


30 


IQ Tau 


5 


3 


5 


10 


20 


(iO 




30 


IRAS 04016+2610 


5 


15 


5 


10 




(iO 


120 


30 


IRAS 04169+2702 




15 


5 


10 


20 


(iO 




30 


IRAS 04248+2612 




15 


26 


10 


20 


(iO 




30 


IRAS 04302+2247 




15 


26 


10 


20 


(iO 




MO 


IRAS 04325+2402 




15 


26 


10 


20 


(iO 




MO 


IRAS 04361+2547 




15 


35 






(iO 


120 


30 


IRAS 04365+2535 




15 


5 


10 




(iO 


120 


30 


IRAS 04368+2557 




100 


100 


10 




(iO 


120 


30 


L1551 IRS5 




15 


15 






(iO 


120 


30 


LkCa 15 


5 


3 


13 






(iO 


120 


30 


RY Tau 


5 


3 


5 


10 


20 


(iO 




30 


UY Aur 


5 


3 


5 






(iO 


120 


30 



Note. — The apertures quoted are in arcseconds. When PSF photometry is done, we took the aperture to be slightly larger than the FWHM 
of the PSF. The apertures for the UBVRIJHKLM data from Kcnyon & Hartmann (1995) are estimates. The 3 aperture for the 2MASS data is 
slightly larger than the 2"pixcl size of 2MASS data. The JHK apertures for DG Tau, DG Tau B and IRAS 04368+25 57 a re those used to carry 
out aperture photometry on 2MASS data. An aperture of 5"was used for the IRAC data from Hartmann ct al. (2005) and Luhman ct al ] 120061 ) 
(the latter used PSF photometry, and 5 is likely an upper limit on the source sizes). The apertures for the IRAC data arc those used for aperture 
photometry by Luhman ct al. (2006) or by the authors of this paper. The MIPS PSF full width half maxima are 6 and 18"for MIPS 24 and 
70 /im respectively. Therefore, we used apertures of 10 and 20 respectively. The apertures on the IRAS telescope were rectangular, therefore we 
are only able to use estimates of an effective aperture radius. Wc used 60 for the 12 and 25 (J.m data, and 120 for the 60 and 100 fj,m data, Finally, 
lAndrews & WilUami j2005l) used a 30 aperture to measure the SHARC II 350 fira sub-mm fluxes, but we do not have information concerning the 
SCUBA and CSO data. Since the FWHM of all these instruments is comparable, wc assumed a 30" aperture for all sub-mm fluxes. 



